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Uber die Restspannungen bei der elastisch-plastischen Biegung 
des krummen Stabes. 


Von W. Swida. 


1, Einleitung. In einer fritheren Arbeit iiber die elastisch-plastische Biegung eines gekriimm- 
ten Stabes wurde lediglich der urspriingliche Belastungsvorgang betrachtet. Die bei der Ent- 
lastung entstehenden Restspannungen und deren Einflu8 auf die gesamte Spannungsverteilung 
bei wiederholten Belastungen wurden dabei nicht untersucht. Dies ist der Gegenstand der vor- 
liegenden Arbeit. 


2. Spannungszustinde. a) Grundannahmen. Wir betrachten einen gekriimmten Stab 
(Abb. 1), der durch die in der Ebene seiner geometrischen Achse wirkenden Momente M bean- 
sprucht wird. Es wird angenommen, daf in derselben Ebene die Symmetrieachse eines beliebigen 
Querschnittes liegt und der Stab vor dem Anbringen der Momente in spannungsfreiem Zu- 
stand war. 


4 € 
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Abb. 1. Abb. 2. Die idealisierte Spannungs-Dehnungslinie. 


Wenn sich unter der Wirkung der Momente im Bereich der groBten Spannungen die plasti- 
schen Zonen gebildet haben, so entstehen im Stab bei der Entlastung, d. h. bei der Verringerung 
der Momente,-Restspannungen. 

Zur Vereinfachung der Berechnung wird ein idealisiertes Spannungs-Dehnungs-Diagramm 
(Abb. 2) angenommen, bei dem die Entlastungslinie O’ A’ parallel zur urspriinglichen Belastungs- 
linie 0 A verlauft. Die Koordinaten o, und &, eines beliebigen Punktes der Geraden 0’ A’ sind die 
Betrage; um die sich die Spannung und Dehnung vermindern. Da auf diese Weise der Zusammen- 
hang zwischen oy, und ¢, durch das Hookesche Gesetz ausgedriickt wird, kénnen deren Werte 
nach den Formeln fiir die elastische Biegung ermittelt werden. 

Weiterhin wollen wir zwei Falle unterscheiden: 

1. Vor der Entlastung bestand die plastische Zone nur an der konkaven Oberflache. 

2. Die plastische Zone bestand sowohl] an der konkaven, als auch an der konvexen Oberflache. 


b) Erster Fall. Wie frither gezeigt wurde?, entsteht bei der Belastung in einer beliebigen 
Faser der elastischen Zone im Abstand y von der Nullachse z die Spannung 
M fl 
Tar ; 


(1) 


o= 

F, 

S,+ 24 (ku—7) 
Ue fal 

Dabei ist S, das statische Moment der Querschnittflache der elastischen Zone in bezug auf die 

Nullinie z, F, die Querschnittsflache der plastischen Zone, r der Krimmungshalbmesser der 


1 W. Swida, Ing.-Arch. 16 (1948) S. 357. 
2 Ing.-Arch. 16 (1948) S. 365. 
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neutralen Flache, y, der Abstand von der neutralen Achse bis zur plastischen Zone und k, die 
Koordinate des Schwerpunktes der plastischen Zone (Abb. 3). 


Die bei der Entlastung verschwindende Spannung in derselben Faser 
ergibt sich auf Grund der oben angefiithrten Annahmen zu 


Ay 

; M ¥y . 
oO aly 2) 

‘ S Tot yy ( ) 


Die Formel (2) ist der bekannte Ausdruck fiir die Spannungen bei — 
der elastischen Biegung. Hierbei ist 2, die Nullinie bei der elastischen — 
Biegung, S das statische Moment des ganzen Querschnittes in bezug auf 
die Nullachse z,, r» der Kriimmungshalbmesser der neutralen Flache 
bei der elastischen Biegung, y, der Abstand der zu betrachtenden Faser 
von der z-Achse. Die Ermittlung von r, und S kann auf bekanntem 
Wege durchgefiihrt werden. Bei vollstandiger Entlastung ist fir M in 
(1) und (2) derselbe Wert einzusetzen. 


4 § : a - ) 
j Die Restspannungen in der elastischen Zone nach vollstandiger Ent- 
Abb. 3. lastung werden: ; 
M if M _ Ye 
Pos Fuy1 TY, Sty Yo : (3) 
Se + ——— (ku—r) 
Tiereyal 

In der plastischen Zone ergibt sich die Restspannung zu 
M Yu 
Oyo =—G 5 4) 

i z S m+» (4) 


wobei o, die FlieBgrenze ist. 

Bei Verwendung der Formeln (3) und (4) sind die Vorzeichen der Koordinaten y und y, zu 
berticksichtigen. 

Die auf diese Weise gefundenen Restspannungen missen die folgenden Gleichgewichts- 
bedingungen befriedigen: 
fooar =o, [ooyaF = 0, (5) 
(F) (F) 
wobei F die gesamte Querschnittsflache ist. Bei der Ermittlung der Spannungen nach (3) und (4) 
kann man die Gleichungen (5) zur Kontrolle beniitzen. 


Fiir einen rechteckigen Querschnitt (Abb. 4) wird 
b 3 b 
S=Fiyiny).  raenyy Faller) =F [at (ry). 


Die Abstande y, und y, lassen sich nach den Formeln bestimmen 


v—t v-+t 


Maa 9 ? Nines 9 Sf P (6) 


Hierbei sind v und t aus den Gleichungen zu errechnen: 
v[A+b(c—v)?] 
PESYOmey (7) 


c v bhv ' 
In = [I Hela (8) 


i 


Os 


ist. Zur Ermittlung der Restspannungen sind die Gleichungen (8), (7), (6) / 
und (3) der Reihe nach zu benitzen. 


Abb. 4. 


c) Zweiter Fall. Fir diesen Fall lassen sich die Restspannungen in der elastischen Zone 
folgendermafen ausdriicken: 


Oy == 


M ¥ ADO Ys 
M1 


. (9) 
S20) r—¥ [Pu ku—r) + Fo (ko +r)] oe 5 Coty) ay 
; 
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_wobei F, bzw. Fy die Querschnittsflachen der unteren bzw. der oberen plastischen Zone (Abb. 5) 
und k, baw. ky die Koordinaten ihrer Schwerpunkte sind. Die anderen Bezeichnungen sind die- 
selben wie im ersten Fall. 


Fir die plastische Zone werden die Rest- 


spannungen 
M Ny 
09 = + Os Sa ree yas (10) 
Das positive Vorzeichen bezieht sich auf die 
obere, das negative Vorzeichen auf die untere 
Zone. 
a ‘ Fur den rechteckigen Querschnitt (Abb. 6) ist 
v 
b 
Se = Se ys) , 
b } 
1 f : r Fy (ku—r) a Q [a?—(r—y,)?] ’ Ser Alen ay 
Abb. 5. Fy (ky ak r) =, 2 [ c?—(r al vg) y Abb. 6. 
Die Abstande y,, y, und r bestimmen sich nach den Formeln 
v—t é vt 2, 
ie ee ee a 
wobei v und t aus den Gleichungen 
eee o= Pop. (12), (13) 
Edy tuy=0, eS 


zu ermitteln sind. Dabei ist 
2 =) lu 


: eto: : (a ee bos ar 3(a+c) ; 
Zur Ermittlung von oy sind (14), (13), (12), (11) und (9) der Reihe nach zu beniitzen. Die 
Gleichgewichtsbedingungen (5) miissen dabei auch befriedigt werden. 


d) Grenzfall. Der Grenzfall, bei dem sich die Plastizierung vor der Entlastung auf den 
ganzen Querschnitt erstreckt, erfordert eine besondere Untersuchung. Bestimmt man auch in 
diesem Fall die Restspannungen nach der Formel (10) und betrachtet r, und S als Kriimmungs- 

halbmesser der neutralen Flache und 

als statisches Moment des gesamten 
Querschnittes, so erhalt man die in : 
Abb. 7 dargestellte Restspannungslinie. \ 
Dann wirden aber die Spannungen im — 
Bereich m n m’ n’ die FlieBgrenze o, C ANS N 
ubersteigen, was bei der angenomme- 
nen Spannungs-Dehnungs-Linie (Abb. 2) | 
nicht moglich ist. Es bleibt daher nur ee | a 
die Annahme ibrig, daB die Zone mn 

Abb. 7. m’n’ von der Hohe y nicht an der Abb. 8. 

Arbeit bei der Entlastung teilnimmt, 
‘d.h. daB sie die bei der Belastung erreichte GréBe ihrer Spannungen o; als konstante GréBe 
beibehalt. Als obere Grenze der Zone mn mn’ muf die Gerade mn betrachtet werden, die den 
Querschnitt in zwei flachenmaBig gleichgroBe Teile zerlegt. Die untere Grenze m’ n’ kann als . 
neutrale Achse eines Querschnittes angesehen werden, der an Stelle des Bereiches m nm’ n’ 
einen Ausschnitt hat. 

Fir den rechteckigen Querschnitt (Abb. 8) lat sich dann die untere Grenze nach der Glei- 


chung bestimmen: 


apes F c—a—yn c—3a+2r, (15) 
TS ee . i a ’ 
dF fo du ‘du yh ean 
U ieee Hf ora si Ear ag 
(F) @ g 


wobei 9 der Kriimmungshalbmesser der Stabachse ist. Nach Ermittlung von ry aus (15) kann 
man auf Grund von (10) die Restspannungen berechnen. 
gts 


\ 


80 Swida: Uber die Restspannungen bei der elastisch-plastischen Biegung. Ingenieur-Archiv _ 
ca nA rah ia a is Na ee eA EL Sea ves EAL I 


Auf diese Weise mu bei der Entlastung eines krummen Stabes eine innere plastische 
Zone entstehen, wenn das Moment bei der Belastung gleich oder nahe dem Grenzwert ist. 

Bei der Berechnung der Spannungen nach Formel (10) ist darauf zu achten, da das Moment 
M, welches bei der Belastung vom ganzen Querschnitt aufgenommen wurde, bei der Entlastung 
in seiner ganzen GroBe von jenem Querschnittsteil zurickgewonnen wird, der in Abb. 7 und 8 
schraffiert ist. Infolgedessen miissen sich die Restspannungen in den einzelnen Punkten des 
Querschnittes etwas von den Spannungen unterscheiden, die man auf Grund des Diagramms 
der Abb. 7 erhalten wiirde. Das tatsachliche Spannungsdiagramm fir ein Zahlenbeispiel ist 
weiter unten in Abb. 12 dargestellt. 


e) Exakte Werte fiir die Restspannungen. Man kann die genaueren Werte fir die 
Restspannungen in einem krummen Stab mit rechteckigem Querschnitt und kreisformiger Achse 
leicht ermitteln, wenn man die friiher abgeleiteten exakten Formeln! fiir die plastische Biegung 
und die bekannten Ausdriicke der Elastizitatstheorie? fir die elastische Biegung beniitzt. Hierbei 


kénnen nicht nur die tangentialen Restspannungen o;, sondern auch die radialen o; bestimmt 


werden. Wenn die plastische Zone vor der Entlastung nur an der konkaven Oberflache bestand, 
so ist in der elastischen Zone 


Oes ees 2 p2 b 2 Ae: Rep aqeen 
O, = Nor 2 r |e a) (1 In -)| + 0 [a(r In . b? In a + rb | 
4 M a” b? b r a 
+ b?(a—o) orl NW (— = In— 1 6? In ; + a? In = + b2—a?), 
oo = — — 5} 5272 eae y+ | a(r? In © + 62 In =) + | 
d Nor iF @ e r 0 (17) / 
+ (0 —a) (b? r?)] 4M (ae ip b 1 62 In r a? in =) i, | 
@ uN Y re a b rie 
wobei 
b b\2 
ey 2 LAS: (Pa ey Bas Cea ead eee 8 
N=b*—9?—20? ln, N’=(b a*)?—4a2b? (In) 
ist. 
In der plastischen Zone wird 
4M 2b? b 
of =o. — ( a In — FB In = + a?ln > + ba?) , (18) 
a 4M / ab b r a 
oe =—o,(1 “) — <3 (4 In +0? In = fatIn =). (19) 


Dabei ist a der Halbmesser der konkaven Oberflache, b der Halbmesser der konvexen Ober- 
flache, r der Abstand vom Kriimmungsmittelpunkt bis zum beliebigen Punkt des Querschnittes 
und @ der Krimmungshalbmesser der Grenzflache zwischen der elastischen und plastischen Zone. 


Wenn sich vor der Entlastung die plastischen Zonen sowohl an der konkaven als auch an der | 


konvexen Oberflaiche gebildet haben, so entstehen nach der Entlastung in der elastischen Zone 
die Restspannungen 


of = — aye {ot? [o (In fe 1 21n > InJ 2) Fa(l | In eer 
+ ao? [rin + of (in © —1)] + 0°[r + 02 (1—In —~)||- (20) 
a ( af In2 ; Pel +a? In & {2 a?) , : 
a — pial en fh tain) bein] “4 
+o ot [o(In ; 1) J #() ln | | gr? [a (in —1) +]! — (21) | 
(28 In es | 62 In ; + @a2In =) ; 


N,'= 0? —o}+20% In, 


1 Ing.-Arch. 16 (1948) S. 361. 


® Siehe z. B. S. Timoshenko, Theory of Elasticity, S. 60, New York und London, 1934, oder A. Féppl und | 


L. Féppl, Drang und Zwang, 1. Bd., S. 305, Miinchen 1924, 


; 
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wobei o und Or die Kriimmungshalbmesser der Grenzflachen zwischen der elastischen und pla- 
stischen Zone sind. Die Spannungen in der unteren plastischen Zone kénnen auch in diesem Fall 
nach (18) und (19) bestimmt werden. 


In der oberen plastischen Zone ergeben sich die Restspannungen zu 


. b 4M a? b? b 
of = or (1—n >) — Wr (— og in + BP In ; ta?in = +52—a?), (22) 
b 4M /a®b? b 
of = —o.In2 — (5 In +8 In— + a?In =), (23) 


f) Restformanderungen. Nun 
ermitteln wir die Restformanderun- 
gen in einem krummen Stab nach 
der Entlastung. Wir betrachten ein 
Element des Stabes (Abb. 9), bei 
dem unter der Wirkung der Momente 
eine plastische Zone an der kon- 
kaven Oberflache auftritt. Wie schon 
friher gezeigt wurde, ist die Zu- 
nahme des Winkels dg unter der 
Wirkung des Momentes M Abb. 9. 


bdo = Mds 


Fay. 
E [S.-+ 222 ke 
Q Faye 
Hier ist 9 der Krimmungshalbmesser der Stabachse. 


A(xg, Yo) 


Die Verringerung des Winkels dy+0ddy bei der Entlastung wird 


Mds vee 
Dann Jat sich die Restanderung des Winkels dq folgendermafien darstellen: 
dodgy = bdy—b.dy = Le it s| (24) 
Gee tera (bs a) 
SBA AL 


Aus Abb. 9 ist zu ersehen, da8 die Drehachsen fiir den Querschnitt pq bei Belastung und Ent- 
lastung verschieden sind. 

Wenn die plastische Zone vor der Entlastung sowohl an der konkaven als auch an der kon- 
vexen Flache vorhanden ist, so wird 


: eee (25) 


| S. Fi os [Fu (ku—r) ai Fy (ko+r)] 
fio hit 


Unter Verwendung eines schon friher angegebenen Verfahrens! erhalten wir folgenden Naherungs- 
wert fir die horizontale Restverschiebung eines beliebigen Punktes A des gekrimmten Stabes 


(Abb. 9): | 
M 1 1 ‘A 

Sti = |bornas + i 

2 | Ee S. any Fu (ky —r) Ss pa 

(%) Fae (26) ‘3 

M 1 1 Nay 

+) / sa | (vey) ds aa 

on) Somes seek [Fu (ku—r) + Fy (ko +r)] Mi 

—¥1 (oa 


Dabei stellt s, die Lange eines Stababschnittes dar, bei dem vor der Entlastung die plastische 
Zone nur an der konkaven Oberflache entstand; s;, ist die Lange eines Stababschnittes, bei wel- ‘ 
chem die plastische Zone vor der Entlastung sowohl an der konkaven als auch an der konvexen on 
Oberflache vorhanden ist. Um den Ausdruck fiir die vertikale Restverschiebung des Quer- 
schnittes A zu erhalten, muB man in (26) (y,—y) durch (x )— x) ersetzen. Der Restdrehwinkel des 
Querschnittes A ergibt sich, wenn man in (26) (y)—y) durch I ersetzt. 


1 Ing.-Arch. 16 (1948) S. 366. 
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3. Beispiele. a) Erstes Beispiel. Ein gekriimmter Stab mit rechteckigem Querschnitt mit 
h—=10 em, b=5 cm, a=10 em, c=20 cm wurde vor der Entlastung durch die Momente M=2tm 
-beansprucht, die die Bildung einer plastischen Zone an der konkaven Oberflache hervorgerufen 
haben. Die FlieBgrenze o; betragt 2400 kg/em?. Auf Grund von (8), (7) und (6) ergibt sich: 

v=9,lom, t=1,9 cm, y,=3,6cm, .y,=5,0 em und r= 14,50 cm. 


Unter Bericksichtigung dessen, daB r, = = 14,43 cm ist, erhalt man das Gesetz der Anderung 
inj 

der Restspannungen: 2 
7260 y 6993 y, 


So 44,5044 14,434, 
y y 


Im Bereich der vor der Entlastung vorhandenen plastischen Zone ist 


6993 y, 
0p = — 2400 Thaw 


In Abb. 10 stellt die Linie ABD den Spannungsverlauf bei der Belastung und die Kurve A’D’ 


den Verlauf der bei der Entlastung verschwindenden Spannungen dar. Die waagerechten Ab- 
stinde zwischen den Kurven sind die Restspannungen. 


Abb. 10, Das Diagramm der Rest- Abb. 11. Das Spannungsdiagramm Abb. 12. Das Spannungsdiagramm fiir das 
spannungen fir M = 2,0 tm. fir M= 2,4 tm. Grenzmoment M = 3,0 tm. 


b) Zweites Beispiel. Derselbe Stab wurde vor der Entlastung durch die Momente 
M=2,4tm beansprucht. Dabei ist eine plastische Zone sowohl an der konkaven als auch an der 
konvexen Oberflache entstanden. Mit Hilfe von (14), (13), (12) und (11) ergibt sich 


&€ = 31,4 cm, v = 8,08 cm, t = 2,08 cm, 
4, = 3,0 ¢m, Vo = 908 cm, r= 14,66 cm. 


Im Bereich der vor der Entlastung vorhandenen elastischen Zone ent- 
stehen die Restspannungen (9) ; 

9340 y 8392 y, 

14,66 14,43-+y, 

Fir die elastischen Zonen ergibt sich nach (10) 


AD 8392 y, 
09 = + 2400 aaa 


Onwia= 


die Linie ABD und der bei der Entlastung verschwindenden Span- 
nungen durch die Kurve A’D’ dargestellt. 
Fur das Grenzmoment ergibt sich nach (15) und (10) 


10510 y 
eS fall} 1) v 
Abb. 13. Die Verteilung der radialen ro = 14,41 cm, O ect en US 14,414 ; 
Restspannungen unter der Wirkung « ‘ Yo 
des Grenzmomentes. Dabei ist bh? ; 
M, = wen 3 tm. 


In Abb. 12 ist das Diagramm von oy fiir den Grenzfall gezeigt. Das Diagramm der radialen | 


Restspannungen of bei der Wirkung des Grenzmomentes ist in Abb. 13 dargestellt. Sie wurden 
nach (19) und (23) errechnet. Man kann sich leicht iberzeugen, da® fiir alle behandelten Beispiele 
die Gleichgewichtsbedingungen (5) befriedigt werden. 


In Abb. 11 ist der Verlauf der Spannungen bei der Belastung durch | 


/ 
| 


i 
4 


at > ; p . ft 
Taks ; 
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4. Der Einflu8 der wiederholten Belastungen. Jetzt nehmen wir an, dai der gekrimmte Stab, 
in dem infolge der urspriinglichen Belastung und der nachfolgenden Entlastung die Restspan- 
nungen entstanden sind, wiederum durch die Momente von gleichem Richtungssinn wie die ur- 
springlichen beansprucht wird. 

Bei Betrachtung eines beliebigen Spannungsdiagramms (Abb. 10, 11, oder 12) ist leicht zu 
erkennen, da infolge des Vorhandenseins der Restspannungen der Vorgang der Formanderung 
bei dieser wiederholten Belastung so lange sich vollkommen elastisch vollziehen wird, bis die 
Biegungsmomente die bei der urspriinglichen Belastung vorhanden gewesene Grée erreichen. 
Bei einer weiteren Zunahme der Momente muB die Plastizierung gleichzeitig in allen Fasern 
entstehen, die bei der urspriinglichen Belastung in der plastischen Zone lagen. Wenn man auf 
diese Weise durch eine voriibergehende Belastung und nachfolgende Entlastung Restspannungen 
in einem gekrimmten Stab hervorruft, so kann man bei wiederholter Belastung von gleicher 
Richtung wie die urspriingliche die GréBe des Biegungsmomentes, das der Stab ohne Auftreten 
plastischer Formanderungen aufnehmen kann, und bei dessen Wirkung das Material sich noch 
-vollkommen elastisch verformt, wesentlich erhéhen. Nach ausgefiihrten Berechnungen betragt 
diese Erhéhung des ,,elastischen’* Grenzmomentes 25 v. H. fiir das erste Beispicl, 50 v. H. fiir das 
zweite Beispiel und 87 v. H. fir den Fall der Belastung durch das Grenzmoment. 

Wenn der Richtungssinn der Momente bei wiederholter Belastung entgegengesetzt ihrem ur- 
spriinglichen Richtungssinn ist, so kénnen (beim Vorhandensein erheblicher Restspannungen) 
schon bei kleinen Werten der Momente plastische Formanderungen in den aufersten Fasern des 
Querschnittes auftreten. So genigt z. B. fiir den in Abb. 12 dargestellten Fall ein Moment von 
0,07 M,, um das Entstehen der plastischen Formanderung an der konkaven Oberflache hervor- 
zurufen. 


(Eingegangen am 21. Marz 1949.) 


Anschrift des Verfassers: Dr.-Ing. Waldemar Swida, (17a) Karlsruhe, Stuttgarter StraBe 21. 
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Beitrage zur Theorie der tragenden Flache*. 
Von N. Scholz. 


1. Einleitung. Die Prandtlsche Theorie der tragenden Linie ersetzt eine Tragflache durch 
einen einzigen Wirbel langs Spannweite, dessen Zirkulation stetig mit der Breitenkoordinate 
der Flache veranderlich ist. Hierbei wird die Strémungsbedingung, daB namlich die resultierende 
Anstrémungsrichtung tangential zur Profilkontur verlauft, nur am Ort der tragenden Linie 
selbst erfiillt. Fir eine Tragflache, deren Tiefe von gleicher GréBenordnung ist wie die Spann- 
weite, ist diese Bedingung unzureichend und ein Ubergang zur tragenden Flache notwendig. Die 
exakte Theorie der tragenden Flache erfullt die Strémungsbedingung in der ganzen durch den — 
FlachenumriB begrenzten Ebene. Sie ist fiir das ebene Problem, d. h. unendliche Spannweite, 
von W. Birnbaum} entwickelt und von H. Glauert? in mathematischer Hinsicht vervollstandigt 
- worden. H. Blenk? versuchte, auch das raumliche Problem der tragenden Flache zu lésen, indem 
er zu den Birnbaumschen Normalfallen fir die Zirkulationsverteilung langs Tiefe noch eine — 
stetige Zirkulationsverteilung langs Spannweite in Form einer Potenzreihe hinzunahm. Die | 
mathematische Durchfihrung bereitet jedoch sehr erhebliche Schwierigkeiten. K. Wieghardt* | 
vereinfachte das Problem zunachst dadurch, da®B er fir die Zirkulation langs Spannweite die | 
elliptische Verteilung annahm und nur die Verteilung langs Tiefe freilie3. Die Stromungsbedin- _ 
gung wurde dabei nur auf der Flachenmittellinie erfiillt. Unter Hinzunahme weiterer Glieder fir — 
Zirkulation langs Spannweite wird dieses Verfahren von J. Ginzel® verallgemeinert. Ferner geht 
ein von K. Wieghardt durchgefihrter Schritt dahin, die Zirkulationsverteilung in Tiefenrichtung 
in einzelne diskrete Wirbellinien zusammenzufassen und die Stromungsbedingung jeweils in der | 
Mitte zwischen zwei Wirbellinien zu erfiillen. V. M. Falkner® ging bei der Behandlung des Pro- 
blems soweit, daB er nicht nur die Zirkulationsverteilung in Tiefenrichtung in diskrete Einzel- 
wirbel aufteilte, sondern nun auch noch die Verteilung langs Spannweite unstetig machte, indem | 
er sie treppenférmig annahm. Die UnzweckmaBigkeit dieses Verfahrens dirfte aber allein daraus 
hervorgehen, daf} damit das sonst stetige Abwindfeld hinter dem Fligel entsprechend den ent- 
stehenden Hufeisenwirbeln eine Anzahl von Unstetigkeitsstellen erhalt, die die Berechnung des 
Abwindes recht unsicher machen. 

Auf Grund dieser Arbeiten erscheint eine Behandlung des raumlichen Problems der tragenden 
Flache aussichtsreich bei Aufteilung der stetigen Wirbelbelegung langs Flachentiefe in diskrete 
Kinzelwirbel, aber unter Beibehaltung einer stetigen Zirkulationsverteilung langs Spannweite, 
die in Form einer Fourierschen Reihe mit freibleibenden Koeffizienten angesetzt wird. Dabei 
soll zunachst im ebenen Fall durch Vergleich mit den exakten Birnbaumschen Lésungen nach- 
gewiesen werden, unter welchen Bedingungen durch die Zusammenfassung der stetigen Wirbel- 
belegung in diskrete Einzelwirbel eine brauchbare Ersatzmethode geschaffen wird. Zumindest 
mu hierbei fiir Gesamtauftrieb und Gesamtmoment der Flache eine Ubereinstimmung gefordert 
werden, die, soweit sie nicht exakt ist, doch schon bei geringer Anzahl von Streifen wenigstens 
gute Naherungen liefert. Ferner ist am Grenziibergang vom raumlichen zum ebenen Fall zu unter-- 
suchen, wie viele Glieder der Fourterschen Reihe fiir die Zirkulationsverteilung langs Spannweite 
notwendig sind, um den richtigen Gesamtauftrieb der Flache zu erhalten. 

Der groBte Teil der fir jeden Einzelfall zu leistenden Arbeit besteht in der Ermittlung der 
Koeffizienten des Gleichungssystems fiir die Zirkulationsverteilung. Diese Koeffizienten stellen 
physikalisch den Abwind in gewissen fest gewdhlten Aufpunkten dar, der ein fiir allemal be- 
rechnet werden kann. Eine wesentliche Abkiirzung der Rechenarbeit fiir den Einzelfall ist dadurch 
zu erreichen, daB dieses Abwindfeld in Form von ,,Abwindtafeln‘‘ universell berechnet wird. 


* Gekiirzte Wiedergabe des theoretischen Teils der Dissertation Braunschweig 1949. (Hauptberichter: 
Prof. Dr. H. Schlichting; Mitberichter: Prof. Dr. Ing. C. Pfleiderer). Der experimentelle Teil befindet sich in 
Forschg. Ing.-Wes. 16 (1949/50) S. 85. 

1 W. Birnbaum, Z, angew. Math. Mech. 3 (1923) S. 290. 

2 H. Glauert, Die Grundlagen der Tragfliigel- und Luftschraubentheorie; deutsch von H. Holl. Berlin 1926. 

3 H, Blenk, Z, angew. Math. Mech. 5 (1925) S. 36. 

4 K. Wieghardt, Z. angew. Math. Mech. 19 (1939) S. 257. 

5 TI, Ginzel, Jahrbuch der Deutschen Luftfahrtforschung 1940, S. I 238. 

° V. M. Falkner, A.R.C. Techn. Rep. a. Mem. Nr. 1910 (1943). 
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2. Die Birnbaumsche Theorie der tragenden Wirbelfliche. Zur Behandlung des ebenen Pro- 
blems der Tragfliigeltheorie hat W. Birnbaum die Vorstellung eingefthrt, die Flache durch eine 
kontinuierliche Wirbelbelegung langs Profilsehne zu ersetzen und kommt somit zur Vorstellung 
der tragenden Wirbelflache. Birnbaum behandelt drei Normalfalle (Abb. 1), fiir die sich Wirbel- 
belegung und Abwind leicht angeben lassen. Durch H. Glauert wurde die mathematische Formu- 
lierung der Theorie mit Hilfe der Methode der Fourierschen Reihen in eine iibersichtliche Form 
gebracht, indem er fiir die Wirbelbelegung langs Flachentiefe die Reihe 


v (9) =V (Cy ctg-S + C, sing +C, sin2p+---) . | (2, 1) 


Wirbeldichte nach der Stre/tenmethode 
bei n=1 2 3 4 Streifen 


+ x 


“@ 
VSB 


. w(x) 


w(a) 


ole) =CoVetg $ = 2azy Vetg $ ng) =G Vsin p =2BVsin p 
L=G FV =10gV T= FVl=Fevl %2(9)-&2 Vsin(2y)- 3-6 Vsin(2y) 
4 2 G y OF: A ie ies CG af rene 
wiV$-£ | v4 fig -2 Pea] w VE beg ael 
lq ~2 Miley (m--F &% Ca = 18 in --#, Cy =9 om - 226 
Erster Norma/fal! Zweiter Normalfall Dritter Norma/fall 


Abb. 1. Die Birnbaumschen Normalfalle. 


ansetzt. Dabei bedeutet V die Anstromungsgeschwindigkeit und y die Zirkulationsdichte. Mit 
dem Koordinatenursprung in der Flachenvorderkante ist der Winkel p definiert durch 


cosp=1+ 75 


und lauft von y=0, x=0 (Vorderkante) bis p=2, x= —I (Hinterkante) (positive x-Achse ent- 
gegen Stromungsrichtung). Der erste Term der Reihe (2, 1) stellt die Wirbelbelegung der ebenen 
Platte dar (Erste Normalverteilung), wahrend die folgende Sinusreihe die Wirbelbelegung in- 
folge beliebig gearteter Wélbung der Flache wiedergibt. Durch Integration von (2, 1) lassen sich 
Auftrieb und Moment einer Wirbelbelegung angeben. Unter Verwendung der Kutta-Joukowsky- 
schen Beziehung A=o VJ" fir den Auftrieb pro Breiteneinheit, wobei die Gesamtzirkulation 


0 
T= f y(x)dx ist, erhalten wir Auftriebs- und Momentenbeiwert! aus (2,1) zu 


x=—l 
A il M rm 1 
Dears oi =2(Cy+5G), Gar On areal 4 (C4 Cn C,). (2, 2) 
2 2 


1 Moment bezogen auf Vorderkante. 
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Diese Ausdriicke enthalten nur die drei ersten Koeffizienten der Reihe (2, 1). Die ibrigen Koef- 
fizienten haben somit keinen Einflu8 auf Auftrieb und Moment. Der induzierte Abwind einer 
Wirbelbelegung ergibt sich nach dem Biot- Oe een Gesetz zu (w>0 Abwind) » 


nije COLE, (2, 3) 


2G 6g 


: FS) 
Die zugehérige Kontur y(x) ist aus der Bedingung zu ermitteln, daB auf ihr die Normalkompo- 
nente der resultierenden Geschwindigkeit verschwindet, so daB der induzierte Abwind 


w(x) =Va(x)= vais . (2,4) 


dx 
wird, wobei «(x) den Neigungswinkel der Kontur gegen die Grundgeschwindigkeit V darstellt. 
Die Gleichung der Kontur erhalt man somit durch Integration von (2, 4). Die wichtigen Werte, 
die sich fir die drei: Birnbaumschen Normalverteilungen ergeben, sind in Abb. 1 zusammen- 


gestellt. 


3. 3 1/4-Theorem und Streifenmethode. Die Ubertragung der Birnbaumschen Theorie der 
tragenden Flache auf den raumlichen Fall bringt vor allem erhebliche mathematische 
Schwierigkeiten mit sich. Es soll deshalb ein vereinfachtes Modell der tragenden Flache ab- 
geleitet werden. Wir betrachten hierzu den weitaus wichtigsten Fall der ebenen Platte und 
ersetzen die zugehorige Wirbelbelegung durch einen tragenden Einzelwirbel+. Dann kénnen wir 


erstens aus der Bedingung gleicher Gesamtzirkulation die Zirkulation des Einzelwirbels an- _ , 


geben zu 


T= [ y(x)dx=na,V1 
; x=—l 
und zweitens aus der Bedingung gleichen Gesamtmomentes die Lage des Einzelwirbels er- 
mitteln: 


1 l 
XV == f xy(a)dz= Re 
x=—l 


Dies ist der vordere Neutralpunkt, fiir den das Moment unabhangig von der Gesamt- 
zirkulation und im Falle der ebenen Platte gleich Null wird. Die Abwindbedingung (2, 4), 
die fiir die stetige Wirbelbelegung langs der ganzen Kontur erfillt ist, kann nunmehr nur noch 
in einem Punkt der Kontur erfullt werden, dem hinteren Neutralpunkt? x4. Die Gleich- 
setzung des Abwindes 


bye 1 
4 al 27 (Ry ee 
4 
ergibt fiir diesen 
3 
x4 Mees 


Wir haben damit gefunden: Wird fiir einen Einzelwirbel in 1/, der Flachentiefe die Abwind- — 


bedingung (2, 4) es */, der Flachentiefe erfullt, so liefert dieser Ersatzeinzelwirbel gleichen 
Gesamtauftrieb und gleiches Gesamtmoment wie die angestellte ebene Platte. 

Die Bedeutung des 31/4-Punktes geht weiter, als es nach dem vorangehenden zunichst er- 
scheint. Wenden wir das 31/4-Theorem auch auf den Fall des Kreisbogenprofiles an, so liefert die 
Abwindgleichung 


iby 
(X24) = a ax ae 
und somit 
r= vi. 


Dyes? stimmt ebenfalls mit dem exakten Birnbaumschen Ergebnis tiberein (vgl. Abb. 1). Wir 
kénnen also ein Kreis bogenprofil durch eine ebene Platte mit gleicher Gesamtzirkulation 
ersetzen, welche im 31/4-Punkt das Kreisbogenprofil berithrt, d.h. einen Anstellwinkel 


1 Die folgende Uberlegung wurde, soweit aus der uns zuganglichen Literatur hervorgeht, erstmalig von 
E. Pistolesi angegeben. . 
2 Nach H. Multhopp auch ,,Anstellpunkt*‘. 
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%g=B/2=2f/l besitzt. Damit ist der 31/4-Punkt des Kreisbogenprofiles dadurch ausgezeichnet, 
das seine Tangentenrichtung mit der Nullauftriebsrichtung zusammenfallt. Das 31/4- 
Theorem erfaBt also im Sinne der Theorie der tragenden Linie neben der Steigung (ebene Platte) 
auch den KrimmungseinfluB der Stromlinien am Ort der Flache. 


Ferner laBt sich zur Theorie des hinteren Neutralpunktes noch zeigen, daB der Schwerpunkt 
der ortlichen Anstellwinkelempfindlichkeit lings Profiltiefe im 31/4-Punkt liegt. Be- 
trachten wir namlich eine langsangestrémte ebene Platte, die nur an der Stelle x iiber eine 
Strecke dx einen 6rtlichen Anstellwinkel x, besitzt, so erhalten wir fiir den Auftrieb langs Tiefe, 
wenn x von Null bis / wandert, 

dca a— 
dx Dalen 


Die Rechnung, die nach der Glauertschen Theorie des Klappenfligels leicht durchfiihrbar ist, sei 
hier ibergangen. Die erhaltene Auftriebsverteilung stellt die der ebenen Platte unter Vertau- 
_ schung von Vorder- und Hinterkante dar und hat somit ihren Schwerpunkt bei 31/4. 


Um nun auch zumindest naherungsweise Angaben tuber die Tiefenverteilung der Zirkulation 
mit Hilfe des 31/4-Theorem machen zu kénnen, teilen wir die Tragflache in eine Reihe hinter- 
einander liegender Einzelflachen gleicher Tiefe auf 
(Flachenstreifen). Jeden Flachenstreifen ersetzen wir IE 
durch einen Einzelwirbel in 1/, der Streifentiefe. Die 
Zirkulation dieses Einzelwirbels ersetzt dann jeweils 
die Gesamtzirkulation des zugehérigen Flachenstreifens 
(Abb. 2). Im Grenziibergang zu unendlich vielen Streifen 
erhalten wir wieder die stetige Zirkulationsverteilung 

der Wirbelflache. Zur Bestimmung der Streifenzirkulation 
erfillen wir jeweils die Strémungsbedingung (2, 4) im P 
31/4-Punkt des Flachenstreifens, was gleichbedeutend ist $ J J 
mit der Mitte zwischen je’zwei Wirbeln. Hierbeiistswenn' |. 7 TSS TTT oe 


es sich um eine gewolbte Kontur handelt, der értliche i; ——}————_.}—- t 


= 4a, ctg 


aneteliwinkel jm 31/4.Dunkt ess otreifens fiir jeden i 8 oe en ¢ 
einzelnen Streifen verschieden. Bei Aufteilung in n Strei- Pe EGE Siero AOS eas RS 
fen haben wir zur Ermittelung der n Teilzirkulationen [, 


2 . 2 x 5 A Abb. 2. Aufteilung der Tragflache nach der 
der Flachenstreifen die Stroémungsbedingungen in n Auf- Siedipeaerhoaal 


punkten zu erfiillen, wobei ein lineares Gleichungs- 
system von n Gleichungen mit n Unbekannten zu lésen ist. Der Gesamtauftrieb ergibt sich: 
nach der Streifenmethode aus der Summe der Teilzirkulationen zu 


Ca Str = = (C, +1,+23 4 ood a). (3, 1) 


Fur das Gesamtmoment um die Vorderkante erhalter wir bei Aufteilung in n Streifen 


—] 
Cm sx = = LT £57 +92, +--+ (4n—3) Ir). (3,.2) 


-Beziiglich der Tiefenverteilung der Zirkulation liefert die Streifenmethode die Wirbelstarke der 
-n. Einzelwirbel, die der Gesamtzirkulation der stetigen Wirbelbelegung in dem zugehdérigen 
Streifenabschnitt entsprechen. Man kann somit eine mittlere Wirbeldichte ftir jeden Streifen 


annehmen von 
Ty 


ies Ue = ate Baw 
Vp, teae y= 1, 2, n ( ) 


Der Schwerpunkt dieser Wirbelbelegung soll ferner jeweils in 1/, der Streifentiefe liegen, so daB 
wir naherungsweise y, als die Wirbeldichte im 1/4-Punkt des »-ten Flachenstreifens ansehen 
_ kénnen. 


4, Anwendung der Streifenmethode im ebenen Fall. Es soll nun zunachst gezeigt werden, 
inwieweit die Streifenmethode im ebenen Fall bei verschiedener Anzahl von Flachenstreifen mit 
den exakten Lésungen nach Birnbaum tbereinstimmt. Dabei ist es zweckmaBig, zur Verein- 
fachung des Gleichungssystems eine dimensionslose Zirkulation zu bilden von der Form 


eas Ie aay ey) 
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Der Faktor 4 ist hierbei eine dimensionslose Bestimmungsgré8e der Profilkontur und wird zweck- 
maBig so gewahlt, daB die dimensionslose Gesamtzirkulation Gges im exakten Falle gleich eins 
wird. Fiir die ebene Platte miissen wir setzen 
Hi Mee 

Dann ist namlich der Gesamtauftrieb nach (3, 1) 

Ca Str = 270g (CG, + G+ --+ Gn) 
und damit das Verhaltnis des Gesamtauftriebes nach der Streifenmethode zu demjenigen der 
exakten Birnbaumschen Lésung 


8 BE (Cg Gea t= Gye (4, 2) 


Ca ex 


Fir das Gesamtmoment um die Vorderkante erhalten wir nach (3, 2) 
1 
Cm Sx = — + og — [Gy +56, + +++ +(4n—3) Gn], 
so daf hier fiir das Gesamtmoment wird 


Cm str = —[6,+56,+ nae + (4n—3)G,]. (4, 3) 


Cm ex 
Wir wollen nun die Gleichungssysteme fiir die ebene Platte bei Aufteilung in verschiedene 
Anzahl von Flachenstreifen aufschreiben. Hierbei stellt jede Gleichung die Stromungsbedingung 
(2, 4) im 31/4-Punkt des jeweiligen Streifens dar. Sie lauten fir n=1 bis 4: 


n=1: Streifenl: : G,=1 
Streifen 1: 26,—26,=1; G= 4 

Pare Seiten, | dig Eis oo Nene 
Streifen 1: 36,—36,— 6,=1; 6,= 7 

n=3: } Streifen 2: G,+36,—36,=1; G,=+ 
Streifen 3: SG AEG et. 3 Coal AC herr (4, 4) 
Streifen 1:  46,—46,—<6,—-+6,=1; 6,== 
Streifen 2: £6,+46,—46,-46,=1; 6,=3 


Streifen 3: 2G6,+26,+46,—46,=1; G,=2 
Streifen 4: 26,+26,+56,;+46,=1; G,=2 


Es ist also jeweils die Gesamtzirkulation nach der Streifenmethode 
6,+6,+---+6,=1 firn=1,2,3,4. (4, 5) 
Dies bedeutet nach (4, 2) exakte Ubereinstimmung zwischen Streifenmethoden und Wirbelflache 


beziiglich des Gesamtauftriebes. AuBerdem erhalten wir, wie sich leicht bestatigen 14Bt, fir alle 
vier Gleichungssysteme (4, 4) 


G, +56, + +.-+(4n—3)G,=1 ‘far n=1,2,3,4. (4, 6) 


Nach (4, 3) bedeutet dies exakte Ubereinstimmung des Gesamtmomentes nach der Streifen- 
methode mit dem Ergebnis der Winbelflache. 


Somit ist fir Gesamtauftrieb und Gesamtmoment gezeigt, daB die Streifenmethode in Ver- 


bindung mit dem 31/4-Theorem fiir die ebene Platte die exakten Werte der Wirbelflache liefert _ 
und zwar bei Aufteilung des Fligels bis zu vier Streifen. Es ist zu vermuten, daB diese Uber-_ 


einstimmung fiir beliebige n gréfSer als vier gilt, ein allgemein gultiger Beweis konnte jedoch 
nicht gefunden werden. Eine Gegeniberstellung der Ergebnisse der Streifenmethode mit den 
Lésungen der Wirbelflache sowohl fir Gesamtauftrieb und Moment als auch fir die einzelnen 
Streifenzirkulationen ist in Zahlentafel 1 gegeben. Dabei sind die exakten Teilzirkulationen der 


einzelnen Streifen zu ermitteln aus ly ex = y(x)dx (Streifentiefe) und diese Werte noch in die 
dimensionslose Zirkulation Gy (4, 1) umzurechnen. 
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Zahlentafel 1. Vergleich der Ergebnisse der Birnbaumschen Theorie mit denen der Streifenmethode 
bei verschiedener Streifenzahl. 


Streifenzirkulation G Gesamt- Gesamt- 
Anzahl ie ay ends Zirkulation Moment 
der 1. Streifen 2. Streifen 3. Streifen 4. Streifen Ca Str/Ca ex = 
Streifen a Tham i ilove ‘i SUN a) 1 mie rie 
(al Sogo eaa MS es 2 6. 3 Biel a | athe y» Cm Str 
n A. KL A «Vi A «nxvl e y, A aVIL | dyer avi Cae 
Str ex Leis? ex Str ex Str ex Str ex 
ce a ae ea a eS El SE eA eile He Ge 
Ebene Platte (erste Normalverteilung) A= 
1 1 1 a — st 1 1 ] 
Re 0,750 | 0,818 | 0,250 | 0,182 c= — — ae 1 1 1 
5: 0,625 | 0,692 | 0,250 | 0,216, | 0,125 | 0,092 — a 1 1 1 
4 0,546 | 0,609 | 0,234 | 0,209 | 0,142 | 0,125 | 0,078 | 0,057 il 1 1 
Kreisbogen-Pro fil (zweite Normalverteilung) 7 =/2 
1 ul 1 — wail il 1/2 
72 0,375 0,500 | 0,625 | 0.500 ae ae a — 1 il 7/8 
3 0,208 | 0,294 | 0,416 | 0,412 | 0,376 | 0,294 aS ae 1 1 17/18 
A 0,136 | 0,196 | 0,293 | 0,304 | 0,317 | 0,304 | 0,254 | 0,196 1 i 31/32 
S-Schlagprofil (dritte Normalverteilung) 1 =36/8 
1 1 0 = — — . — — — 1 0 —] 
By 0,562 | 0,424 |—0,313 | —0,424 1/4 0 1/2 
3 0.347 0,370 | 0,139 0 (nd Oi —Oss LO — — 1/9 0 7/9 
4 0,239 -| 0,276 | 0,220 | 0,148 0,079 | —0,148 |—0,317 | —0,276 1/16 0 7/8 


In der gleichen Weise kann nun auch das Kreisbogen- und S-Schlagprofil (zweiter und dritter 
Normalfall) behandelt werden. Beim Kreisbogenprofil bilden wir eine dimensionslose Zirkulation 
nach (4,1) mit A=§/2 und beim S-Schlagprofil mit A=36/8. Dann lauten namlich die Abwind- 
gleichungen in den so definierten dimensionslosen Zirkulationen Gy ganz analog zu denen der 
ebenen Platte (4, 4), wie man sich leicht iiberzeugen kann. Lediglich das inhomogene Glied der 
rechten Seite ist jetzt nicht mehr éins wie bei der ebenen Platte, sondern nimmt den Wert «(x)/A 


an, der im wesentlichen den 6értlichen Anstellwinkel der Kontur im 31//4-Punkt des jeweiligen ~ 


Streifens darstellt. Der Wert «(x) ist aus den Abwindgleichungen der Birnbaumschen Normal- 
verteilungen sofort zu erhalten. In Zahlentafel 1 sind die Lésungen dieser beiden Falle ebenfalls 
enthalten. Bemerkenswert ist vor allem, da sich auch fir das Kreisbogenprofil exakte Uber- 
einstimmung des Gesamtauftriebes nach der Streifenmethode mit der Birnbaumschen Lésung bei 
Aufteilung in ein bis vier Streifen ergibt. Das Gesamtmoment liefert bei drei bis vier Streifen 
bereits geniigende Annaherung an den exakten Wert, so daB eine Streifenzahl von 3 bis 4 als 
ausreichend betrachtet werden kann. 

Wir sehen hieraus, da® die Streifenmethode in Verbindung mit dem 31/4-Theorem Gesamt- 
auftrieb der ebenen Platte und des Kreisbogenprofiles sowie Gesamtmoment der ebenen Platte 
exakt liefert. Fiir das Gesamtmoment des Kreisbogenprofiles sowie des S-Schlagprofiles erhalten 
wir bereits bei einer Streifenzahl von 3 bis 4 recht gute Naherungswerte. Dartiber hinaus nahern 
sich die Teilzirkulationen der einzelnen Streifen den exakten Lésungen fir jede Kontur mit zu- 
nehmender Streifenzahl. Hierbei liegen die Teilzirkulationen der vorderen Wirbel durchweg unter 
den exakten Werten, die der hinteren dariiber. Das gleiche geht aus Abb. | hervor, wo jeweils im 
31/4-Punkt des Flachenstreifens die GréBe der Wirbeldichte nach (3,3) in die Birnbaumsche 
Wirbelbelegung eingetragen ist. Da fiir,die Birnbaumsche Theorie das Uberlagerungsprinzip gilt, 
darf geschlossen werden, daB auch die Streifenmethode fiir eine beliebige Kontur bei Aufteilung 
in wenige Streifen recht gute Naherungswerte liefert. 


5. Die Ubertragung auf das rdumliche Problem. Fiir das raumliche Problem tritt nun eine 
weitere Abhangigkeit des Abwindes und damit auch der Zirkulation von der Breitenkoordinate y 
der Flache hinzu, so daB die grundlegende Strémungsbedingung, die im ebenen Fall durch (2, 4) 


gegeben war, jetzt lautet w(x,¥) = Vox), SA) 
was auf der ganzen Tragflache zu erfillen ist. Der Abwind setzt sich dabei infolge des hinter der 


Flache abgehenden freien Wirbelbandes aus zwei Anteilen zusammen, dem gebundenen und dem 


freien Anteil: w= w+ uy. 
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Diese Anteile lassen sich nach dem Biot-Savartschen Gesetz angeben. Wir berechnen sie fur einen 
geraden tragenden Wirbelfaden, der in der y-Achse liegt (Abb. 3). Seine Zirkulation an der 
Stelle y’ sei [’(y’), dann hat der abgebende Wirhelfaden die Starke (dI'/dy’). Der Abwind in 
einem Punkte (x, y) der Flachenebene betragt, wie aus Abb. 3 abgelesen oe kann, 


+b/2 
{I 4 cos py , 
W,(%,¥) = 4 ii Py.) epee ye » (5, 2a) 
y= —b/2 
+b/2 ee 
1 di’ (1+ cosy) i; 
wy (X5 Y) ery dy’ (y—y) dy (5, 2b) 
Ff y’ =—b/2 
\ re: 
\ mit tgy= (anys 


Wenden wir jetzt auf unsere Tragflache die 
Streifenmethode an, wobei wir uns wegen der Grad- 
linigkeit der gebundenen Wirbelfaden auf den recht- 
eckigen Grundrif§ beschranken missen, so summieren 
sich die Abwindanteile der einzelnen Streifenwirbel. 
Legen wir die y-Achse in die Flachenvorderkante, so 
haben die Lagekoordinaten der Streifenwirbel ["(), die 
jeweils in 1/, der Streifentiefe liegen, die Werte 


(=x) =F 1 (y=1,2,8,...n). 


llbvstststdedidiadatsedediae Vides 


Abb. 3. Zur Aufstellung der Abwindformeln. 


Dann ergibt sich aus der Perna (3,1) die fundamentale Gleichung des Problems 


+b/2 
J n () cos wy ssa eug (1 + cos <0) 


y= —b/2 y= —b/2 


a (x,y) = 


7 


mit tgy = aaah (Di Ue Zed valet) 


Xy— X 


(2) 


mit den gesuchten Funktionen ly: Lys Py); - : . Dy). Diese Integrodifferentialgleichung ist 
zu erfiillen auf den Linien x=konst zwischen —b/2 ee +6/2, welche die 31/4-Linien der einzel- 


nen Flachenstreifen darstellen, also fiir 

(Ge Se 1.2.8.2 a). 
Gleichung (5, 3) kann fiir n=1 (Traglinientheorie) und verschwindende Tiefe 1<b in die Prandtl- 
sche Integrodifferentialgleichung tberfiihrt werden. Hierzu miissen wir noch das vom gebundenen 
Wirbel herrithrende Integral durch partielle Integration etwas umformen. Wie in Abschnitt 8 
noch kurz a serciey wird, erhalt man 


4/2 
ie ue inpdy’ 
Ext A 7 (—<x) dy’ Bir sae (ite 
y= —b/2 


Fir unseren Grenztibergang //b—>0 und n=1 haben wir zu setzen 
(x) —(—x,)=l/2; p=n/2 fir y’>y und y=—n/2 far Y <uyi 


Dann laBt sich das erste Integral ohne weiteres auswerten, und wir erhalten die bekannte Prandtl- 
sche Gleichung 
I(y) 1 ye died 
«(y) a VI(y) 1 aoe, J a dy yy" 
aol 

Der erste Term dieser Gleichung stellt Leech. ee effektiven Anstellwinkel «, dar, whe ist 
im Sinne des 31/4-Theorems der Abwind des gebundenen Wirbels auf der 31/4-Linie der Trag- 
flache bei unendlich langer tragender Linie (ebener Fall). Der zweite Term stellt den Abwind des 
freien Wirbelbandes auf der 1/4-Linie, also'am Ort der tragenden Linie selbst dar. Die beiden 
Terme enthalten gegeniiber dem 31/4-Theorem Abweichungen im entgegengesetzten Sinne, so 
daB sie sich teilweise kompensieren. Hierauf diirften die auch bei kleinen Seitenverhaltnissen 
noch recht guten Ergebnisse der Traglinientheorie zurtickzufiihren sein. 


(5,4) 
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Um eine Lésung der Gleichung (5, 3) zu finden, vereinfachen wir das Problem jetzt noch da- 
durch, daB wir diese nicht in allen Punkten der 31/4-Linien x=konst langs Spannweite erfillen 
sondern nur in einzelnen Punkten. Dann laBt sich namlich die Integralgleichung (5, 3) zerlegen 
in ein lineares Gleichungssystem von endlich vielen Gleichungen. Fiir die Zirkulationsverteilung 
I(y) erhalten wir damit Naherungen, die noch einer genaueren Betrachtung unterzogen werden 

sollen. 
Wir setzen hierzu die Zirkulationsverteilungen der einzelnen Streifen als die folgenden tri- 
gonometrischen Reihen an: 
/ (1) : 
Py) = V bag (a, sin 8 + a; sin3 9 +a,sin59+.--), ; 
(2) (: > 5) 


Dabei bedeutet a, den geometrischen Anstellwinkel an irgendeiner fest gewahlten Stelle der 
Flache, auf den alle Neigungswinkel der Flache bezogen werden. 
Die Winkelkoordinate @ definieren wir in bekannter Weise durch 
| cos } = es, ‘ 
so daB in den Grenzen von + 6/2 bis — b/2 der Winkel # von 0 bis a variiert. Da die Funktion 
I’(y) symmetrisch zur Flachenmittellinie ist, scheiden die geraden Glieder von vornherein aus. 
Fiir den Abwind eines Punktes der Flachenebene machen wir dann entsprechend (5, 5) den Ansatz 


_ e =e) + e2) 4 6) 4... 4 elm) (5,6a) 


oder 


€ = og [ (a; ef + ag ef) + ase +---) + (be + bye) + bee) +--+) 4 -+-], (5, 6b) 


wobei die el), et en .. als noch zu berechnende Abwindfunktionen aufzufassen sind, die sowohl 
von der Breitenkoordinate y als auch von dem relativen Abstand zwischen Wirbel und Aufpunkt 
(x—x,) abhangen. Zur Ermittlung der Koeffizienten a,, a3, a;,... b,, bg, b;.... usw. steht uns die 
Stroémungsbedingung (5, 1) in der Form 

E(x, y) = (x,y) (5,7) 


zur Verfiigung. Wir kénnen diese also bei n Streifenwirbeln und m Koeffizienten der Reihen (5, 5) 
in n-m Punkten der Tragflache erfiillen, woraus sich ein Gleichungssystem von n-m Gleichungen 
ergibt. Aus Symmetriegrinden kénnen wir diese Aufpunkte nur auf einer Halfte der Flache 
wahlen und nehmen hierbei aquidistante Abstande des Winkels 0. Unser zu lésendes Gleichungs- 
system nimmt dann die folgende Form an: 


a, ef) + ag ef) + a, ef) + >> | 
by eb) + by ef) + Bs es) + ++ UA eC (5,8) 
(3), (3) (3) OS Zeiger der Aufpunkte). 

C) €} 5 1 €3 35.1 os 55 1" Braet 
Das Koeffizientenschema der linken Seite dieses Gleichungssystems ist nur abhangig vom Seiten- 
verhaltnis bzw. der Flachenform, nicht dagegen von der 6értlichen Anstellwinkelverteilung (Ver- 
windung und Wolbung). Diese wird nur durch die inhomogenen Glieder «,/xg bestimmt. Als 
Lésung des Gleichungssystems erhalten wir die n-m Koeffizienten der Gleichungen (5, 5), die uns 
die Zirkulationsverteilungen der einzelnen Streifenwirbel langs Spannweite liefern. Die Auf- 
triebs- und Momentenbeiwerte der Flache schreiben sich mit diesen Lésungen, wie folgt [vgl. im 


ebenen Fall Gleichung (3, 1) und (3, 2) |: 


ae ye [(a, +b, +e,4+-+-) sind + (ay+b,+e,+---) sind d+ (5, 9) 
+ (a;+b,+e,+---) sin5 O+-:-] hg » 

Ca ges = 7 a (a,-++5,-+¢,+ 30 -) Ag» (Sy. 10) 

emont = — go P [(a, + 5b, +9e, +--+) sin + (ap +50, +94) sim39+| Os Vy) 


+ (ag +5b,+9c,) sin5 O+---] a, 
‘ EF (a $55, +9 +++) ty: (5,12) 
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Hierzu ist bemerkenswert, da® zu den Gesamtwerten von Auftrieb und Moment nur die ersten © 


Koeffizienten der Reihen (5,5) a,, by, c,,... Beitrage liefern. Diese ersten Glieder stellen die ellip- 
tische Verteilung dar. 


6. Die Grenzfalle der Zirkulationsverteilung lings Spannweite. Wir haben die Zirkulations- 
verteilung langs Spannweite als eine unendliche trigonometrische Reihe (5, 5) angesetzt, deren 
Gliederzahl aber zunichst noch offen gelassen wurde. Wieviele Glieder der Reihe (5, 5) not- 
wendig sind, um bei ertraglichem Rechenaufwand eine ausreichende Genauigkeit zu erzielen, soll 
am Grenztibergang zum ebenen Fall nachgeprift werden. Hierfiir lassen wir —+0 gehen, so daB 
b/l—> © geht. In diesem Falle liefert nur der gebundene Wirbel einen Anteil zum Abwind. Dieser 
betragt auf der 31/4-Linie mit dem Abstand 1/2 vom gebundenen Wirbel nach dem Biot-Savart- 


schen Gesetz 


etn > pea 
&g(y) anv nlV~ 


Da I—0 geht, multiplizieren wir mit J/b und erhalten mit ["(y) nach (5,5): 


9 : ty Hy pate =é (a, sin }-+-a, sin3 Pa, sind + ---). 
Bedenken wir, daB unser Ansatz fiir den Abwind nach (5,6b) lautete =a (a, €, +43 €3+45&;+-+-), 


so folgt sofort aus einem Vergleich 


l : l f l : 
7% > & = sin; 7&3 = sin3 B; Mb = sin5;.... (6, 1) 


Es ist nun das Gleichungssystem (5, 7) zu lésen fir den Fall der ebenen Platte ohne Verwindung, 
also mit %g =p. Damit ist zunachst das inhomogene Glied fiir alle Gleichungen gleich eins. Ferner 
kénnen wir uns wegen //b—>0 auf einen tragenden Wirbel beschranken. Erweitern wir die linken 
Seiten der Gleichungen (5,7) noch mit dem Faktor z1/b und fihren ein 


at — —“# (M21, 3/5 oe) (6, 2) 


Pe en gen TO 
wobei aj, fiir 1/b—+0 endlich bleibt, so erhalten wir das folgende Gleichungssystem: 
1 abe are ; 
ay (x 7 e1p) | a3 (x b eap) | ae (x 5 esp) Soca | 


(p=1, 2, 3, ...m, Zeiger der Aufpunkte) . 


(6, 3) 


Die Klammerausdricke haben wir bereits in (6, 1) ermittelt. Die Aufpunkte wahlen wir auf der 
31/4-Linie in aquidistanten Abstanden des Winkels 3? zwischen #=z/2 und #=0. Dann lauten 
die Gleichungssysteme fir verschiedene Gliederzahl m: 


m=l: b= Fi aj = 1, ay = 1, 
h 
b= >: ay ae a3 ae aha >+ ye, 
Ws 
1 1 1 1 1 
1 
OS ha ab ce of tel + _ 2 
2 *, OI ea i Vea 
V3 V3 
Wider) b=: yet + 0-43 : ig ee Ai, aba, 
Atay rn * * A * i ee 2 ul 
o 6° ) ay } az + 2 Cra, Ge Gane Ions 


Fur den Auftriebsanstieg der Tragflache erhalten wir nun, indem wir (5,9) mit l/b erweitern, 


deg Tay at 
i] 1 
dug 2 Te 2 


und damit fiir das Verhaltnis des Auftriebsanstieges der Naherungslésung zu dem exakten Wert 
27 im ebenen Fall 


1 aloe BB se } 
F ae ginerlne ; 
22% da&go 4 


R32 : 
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Dies ergibt die folgenden Zahlenwerte: 


a ae eS En 
Tt pg Ml 1 nm (2 1 

eI MR en aa 
4 4 \2 4 7 

2m da&ge EW V2 3 V3 
0,786 0,945 0,978 


a eeeeeeeeeeSSSEeee 


Bei drei Sinusgliedern betragt danach der Fehler noch 2°, was als ausreichend angesehen 
werden kann. Das Ergebnis bei m=1 erhalt K. Wieghardt im Grenziibergang zum ebenen Fall, 
da er sich auf das erste Glied, die elliptische Verteilung, beschrankte. Die Zirkulationsverteilung 
langs Spannweite schreibt sich nun mit J’, = Vimag,, in den Koeffizienten ai, 


Ly) 
Do 
Sie ist in Abb. 3 fir m=1, 2, 3 tber der Spannweite aufgetragen. Fiir m=3 ist die Annaherung an 


die konstante Auftriebsverteilung des ebenen Falles bereits recht gut. Dabei mége bedacht wer- 
den, daB hier der ungiinstigste 


= aj sin? + a3 sin3 } + a¥ sind P+... . 


Fall vorliegt, ein Rechteck Ty) 

durch eine endliche Sinusreihe = 
anzunahern. Bei den Flachen ee 
von endlichem Seitenverhalt- 

nis kommt jedoch die Auf- ci 


triebsverteilung der ellipti- - 
schen naher, und man kann 
deshalb erwarten, da in sol- 
chen Fallen die Annaherung 
durch die dreigliedrige Fourier- 
Reihe wesentlich besser ist als 


nach Abb. 4. 


Betrachten wir nun noch 
den Grenzfall eines sehr klei- 
nen Seitenverhaltnisses b—0, 
so rickt die 3//4-Linie sehr 
weit hinter den tragenden 
Wirbel, so daf nur die freien 
Wirbel Beitrage zum Abwind liefern. Fir den Fall der unverwundenen Tragflache mu der Ab- 
wind auf der 31/4-Linie konstant, das heiBt unabhangig von y sein. Fir x—>— ~ besitzt nur die 
elliptische Zirkulationsverteilung einen konstanten Abwind, namlich 


Abb. 4. Annaherung des ebenen Falles durch m Glieder einer Sinusreihe 
fir die Zirkulationsyerteilung langs Spannweite. 


Ca 

ai © 

Die unverwundene Tragflache sehr kleinen Seitenverhaltnisses besitzt also unabhangig von der 
FlachenumriBform eine elliptische Auftriebsverteilung. Aus éu=ag erhalt man den Auftriebs- 
anstieg 


X—>— 0: &y = 2 


dca ss 
dy = 94. 


Er ist gerade halb so groB wie der der Traglinientheorie bei A—>0. * 


A> 0% 


7. Die Zirkulationsyerteilung lings Flaichentiefe. Aus der Streifenmethode erhalten wir eine 
Zirkulationsverteilung langs Tiefe nur insofern, als wir die Gesamtzirkulation der einzelnen 
Flachenstreifen kennen. Am ebenen Fall wurde bereits gezeigt, daB drei bis vier Streifen eine aus- 
reichende Annaherung an die exakten Werte der Wirbelflache bringen. Wir wollen jedoch aus 
_den Ergebnissen der Streifenmethode noch eine stetige Zirkulationsverteilung tber die Tiefe er- 
rechnen, um somit die Druckverteilung tiber der gesamten Tragflache angeben zu kénnen. Hierzu 
setzen wir fiir die Verteilung langs Tiefe die vom ebenen Fall her bekannte Glauertsche Reihe an, 
in deren freie Koeffizienten wir jetzt die Anderung der Zirkulationsverteilung langs Spann- 


* Fiir die Traglinientheorie ist dc,/dag = 20 A/(2+ A). 


7 4 Sy SE ie eg ee a, ee Ph \s 

“i Oats Soe os pt Beas 
x . Nae “ 
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weite legen. Wir erhalten also eine flachenhaft tiber den Flachengrundri8 verteilte Zirkulation in | 


der Form 


y (x.y) = Vtg [(ko: sin O + kgs sin 3.9 + igs sin 9 +--+) etg + | 


| 
| 


+ (k,, sin? + k,, sin3 } + ky, sind @ + ---) sing + | (7, 1) 


+ (kg, sin 0 + keyg sind & + kegs sind? + .--)sin2g+-- ale 
wobei die Winkelkoordinate m die Anderung langs Tiefe und die Winkelkoordinate &@ die Ande- 


rung langs Spannweite liefert. 
Zur Bestimmung der Koeffizienten in (7, 1) benutzen wir die Bedingungen gleichen Auftriebes 
und gleichen Momentes wie bei der Ser oreee fir jeden Flachenschnitt Idy, also 


(24) ev frtanae= (44)... (7, 2a) 
ao 

(2) mov fx rtanaem (8), (7, 2b) 
" x=—l 
0 

or han OU Pee aaa (7,26) 
x=—l 


Dabei entnehmen wir die Werte der rechten Seiten aus den Lésungen der Streifenmethode, wobei | 


je nach der Streifenzahl n dieser Gleichungen (7, 2) benétigt werden. 

Das Verfahren liefert im Grenziibergang zum ebenen Fall dann die exakt richtige Zirkulations- 
verteilung uber die Tiefe, wenn auch die Streifenmethode die richtigen Werte fir Gesamt- 
auftrieb und Gesamtmoment im ebenen Fall liefert (vgl. Zahlentafel 1). In den Fallen, in denen 
die Streifenmethode von der exakten Lésung abweicht, mu} das Ergebnis der Streifenmethode 


mit einem Korrekturfaktor (Mex/Mstr)45. versehen werden, um nach (7,1) einen Ubergang in | 


die exakte Tiefenverteilung des ebenen Falles zu erhalten. 


8. Die Berechnung der Abwindfunktionen. Wir haben bisher noch die Berechnung der Abwind- 
funktionen e(x, y) offengelassen, die wir als Koeffizienten des Gleichungssystems (5, 7) benétigen. 
Diese sind universell, und ihre Berechnung nimmt damit den Hauptteil der fiir einen speziellen 
Fall notwendigen Rechenarbeiten vorweg. Wir fihren sie durch fiir den Fall eines geraden 


tragenden Wirhels, dessen Zirkulation J’(y) nach den ersten drei Gliedern der Sinusreihe (5, 5) 


verteilt ist, also fir w=1, 3, 5 (vgl. Abb. 3), 

Den Abwind in einem Punkte (x, y) der Tragflache hatten wir bereits in (5, 2) aufgestellt. 
Der Abwind der freien Wirbel enthalt fiir Punkte hinter der Tragfliche (x—0) eine Singularitat 
bei y=y’. Um wy bequemer berechnen zu kénnen, spalten wir daher fir x <0 das Integral (5, 2b) 


auf in 
+b/2 +b/2 
1 di’ dy! 1 dI’ 1—cosy 5 


ee dy’. (8, 1). | 


ln diy 9 An dy y—y’ 
=b/2 —b/2 


Jetzt liegt die Singularitat nur noch im ersten Integral, das sich nach H. Glauert analytisch aus- 


werten laBt. Das Abwindintegral der gebundenen Wirbel laBt sich als eine Integration tiber p 


schreiben. Nach Erweiterung mit —x wird mit dy’ = —xdy/cos*yp nach (5, 3c) 
+b/2 
Ws = wos jbnel ‘) cospdy. 
—b/2 

Durch partielle Integration spalten wir dieses auf in 

4/2 4/2 

ne are ; 
er copa {[Po ) sin y| — i dg sinpdy'. 
—b/2 =b/2 


Der erste Term wird Null, da J"(y’) an den Grenzen verschwindet, so daB wir endgiltig fiir dent | 


Abwind der gebundenen Wirbel erhalten 
ee 


We a [3 9  iny dy’. (8, 2) 


y= —b/2 


ned 
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Wir fihren nun mene den dimensionslosen Koordinaten 7 = ae und = ah fir die Zirkulations- 
verteilung I’(y) den Reihenansatz (5, 5) ein 

I(y’) = Vbag(a, sin ® + a; sin3 & + a, sin5 ®”) 
und ebenso fiir den Abwind eines Punktes der Tragflache nach (5, 6b) 

{ E = (4, & + 4&5 + As Es) ag 


und bei Aufteilung in den Beitrag des freien und des gebundenen Wirbels 


es La, (ef, + &g1) + @3(é3 + Egs) + 45(é5 + zs) | Kgs 


a) Freier Wirbel: 

4 = sn(u®) 1 1, (l—cos y)d & 

4 peter sin ? aa Jo tae ) cos 0 —cos ? (u=1, 3, 5), (8, 3) 
z b) Gebundener Wirbel: — e,, = pee cos (uH) sinpd® (4=1y3) 8), (8, 4) 
_ ¢) Gesamtabwind: ¢, = [i ae 5 le J 0 (wd) tg + dé" (8, 5) 
3 9’=0 


mit tg p= cos 0 — cos } 


=e 

Far die Grenzwerte bei €—+0 (am (Elegy 
- Ort des Wirbels) erhalt man fi 
pe ‘ I 
‘ ee snGey C1", ae Re | 
C/E anne ae Ae ere 

e sin (8) L geb. Wirbe: Asal 
eee .(80b) al i 
-und bei —»— co (sehr weit hinter Ta 

dem Wirbel) a1 

; sin (uw): i 

SERS eal tee a |(@, 60) l ee) 
a Ue -20 -15 —-10 —05 05 40 15 20 


Ferner findet man folgende Sym- 
-metriebeziehungen: 


es ry) 
Be oy ay h8) | ell : 


a. (5 bar = 8 i (+6), (8, Tb) 8 


(ey | 
ew) 0, (48). | (7) 


| Damit beschrankt sich die Berech- 
-nung der Abwindfunktionen auf 
| positive Werte von € und wegen der 
| Symmetrie zur Flachenmittelebene 
| auch auf positive Werte von 7. 

Die in den Formeln (8,3), (8, 4) 
-und (8,5) auftretenden Integrale 
| sind von elliptischem Typus, lassen 
} sich aber im allgemeinen nicht auf 
| elliptische Integrale erster und zwei- 
» ter Gattung zuruckfiihren. Von den Abb. 5. Abwindfunktionen des ersten Fouriergliedes der Zirkulationsverteilung. 
B* 


freie Wirbel 


| 
Dy 
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+1 Nullstellen der Integranden werden yu Nullstellen durch das Glied cos w@ verursacht, 
und eine Nullstelle durch den Winkel y, der bei o =9 «zu Null wird. Diese Nullstelle verdient| 
besondere Beachtung, da im Falle £9 an dieser Stelle eine Unstetigkeit mit Vorzeichenwechsel 
auftritt. Fir kleine € und bei ~=1 wurde zur numerischen Integration das Integrationsintervall 
in drei Bereiche aufgeteilt, namlich in 


OSS a Os 1 Ss) 


AN 


rt aa j 
(- S)egs Fir €>4 konnte die Aufteilung ; 
—s 1 bei ® = % unterbleiben. Fiir jedes 
ye as Teilintervall wurde die Quadra-: 
tur mit der Gaufschen Formel! 
geb. Wirbel bei fiinf Soe durch- > 
+ 7 gefihrt. Fir die Berechnung der’ 
Integrale mit w=3,5... laBt . 
sich das w-fache des Arguments | 
des cos zurtickfiihren auf eine 
Summe von Potenzen des cos, 
wobei die héchste Potenz vom. 
Grade mw ist. Hierdurch behalt 
der Integrand in allen Fallen 
seine Charakteristik mit zwei. 
Nullstellen bei, so da die Qua-) 
dratur ganz analog durchgefihrt | 
werden kann. Die Genauigkeit | 
des Endresultates wird allerdings | 
durch die anschlieBende Differenz- 
bildung der Einzelquadraturen 
verringert. Es wurden im allge- 
meinen nur die Abwindanteile ¢, ,, 
(Abwind des gebundenen Wirbels} 
und ¢, (Gesamtabwind) berech- 
net und nur zur Kontrolle auch 
einige Werte von e, (Abwind 
der freien Wirbel). Zur graphi- 
schen Darstellung von e,, und 
freie Wirbe/ é, ist es zweckmafig, das Pro- 
dukt (—&)e, tiber € aufzutragen. | 
da dieses fir —>+0 gegen ein) 
endlichen Wert geht. Der Verlauf 
der berechneten Abwindfunktio- | 
nen ist aus Abb. 5 bis 7 zu er-| 
sehen. Die Zahlenwerte der Ab-. 
windfunktionen befinden sich im) 
Anhang in Zahlentafel 5. Dort | 
ist ebenfalls ein praktisches 
Rechenbeispiel fiir die Verwen-. 
dung der Abwindtafeln aufge- | 
fiihrt (Zahlentafel 6a, b, c¢, d). | 


Abb. 6. Abwindfunktionen des dritten Fouriergliedes der Zirkulationsverteilung. 


9. Anwendungsbeispiele. a) Elliptische Flache. Wir wollen zunachst die Flache mit| 
elliptischem Grundri®B bei gerader 1/4-Linie betrachten, die wir nach dem 31/4-Theorem durch 
einen Einzelwirbel in der 1/4-Linie ersetzen, wobei die Strémungsbedingung auf der 31/4-Linie 
erfillt wird. Setzen wir elliptische Zirkulationsverteilung langs Spannweite voraus, wie dies die 
Traglinientheorie liefert, so haben wir zur Bestimmung des einen freien Koeffizienten a, der 
Zirkulationsverteilung die Abwindbedingung (5,7) in einem Aufpunkt der 31/4-Linie zu a 


fiillen, welche dann lautet 


- Pe en Ay pe a aA Ute er eae a Fy SH cided \ “lie 
‘ PA Rr Pen en tae E e 
ae Rare Yaad ers e 
Sa "a5 Peay 
es } 
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Mit (5, 10) erhalten wir fir den Auftriebsanstieg 


dtq nA al 
dag nr ts (9, 1) 


Den Aufpunkt wahlen wir im Symmetrieschnitt auf der 31/4-Linie der Tragflache, wo sich die 
Abwindfunktion ¢, in die Form tberfiihren labt 


1 ee 1] 3) 


TUS lal ey 


(E ist das elliptische Integral zweiter 
Gattung). Der Abstand (—&) in die- 
sem Punkte ist ftir die elliptische 


Flache 
4 


( E) 31 a tea nA’ 
so da wir aus (9, 1) erhalten 
dea 
dug ee (9, 2) 


22A 


a+ eG) 8 (ar z) 
IG ES 2 


Diese Formel laBt sich vereinfachen, 


wenn wir fiir das elliptische Integral 
die Naherung einfiihren 


2 
1 fa eG é) ee 
E (———., 3) | 
V B11 2. ft] 
Damit erhalten wir die auch von 
H. Schlichting in einem unverdffent- 
lichten Bericht auf andere Weise ge- 
wonnene Formel 
don. 272A 
dw 94 A214 


(9, 3) 


Die Prandtlsche Formel erhalten wir, 
wenn wir fiir den Abwind die fiir 
groBe Seitenverhaltnisse giiltige Nahe- 
rung machen 


Ld 
Ease ons bete, Gas 
Dann ergibt sich namlich aus (9, 1) 
dey 227A 


(9, 4) 


begs 


Abb. 7. Abwindfunktionen des fiinften Fouriergliedes der Zirkulationsverteilung. 


Einen Vergleich dieser Formeln zeigt Abb. 8, wo auch vergleichsweise die Werte fir den Ellipsen- 
fliigel nach K. Krienes! eingetragen sind. Man kann den Auftriebsanstieg des Ellipsenfliigels auch 
statt mit einem mit drei Aufpunkten langs Halbspannweite nach dem 31/4-Theorem ermitteln 
und erhalt dann die ebenfalls in Abb. 8 eingetragenen Werte, die nur wenig von der Kurve 
nach (9,2) abweichen. Dies besagt nichts anderes, als daB auch nach dem 31/4-Theorem der 
Ellipsenfligel naherungsweise eine elliptische Zirkulationsverteilung besitzt. 


1 K. Krienes, Z. angew. Math. Mech. 20 (1940), S. 65. 
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Zahlentafel 2. Koeffizienten der Zirkulationsverteilung | 
lings Spannweite, Auftriebsanstieg und Neutralpunktlage der ebenen Rechteckplatte 
nach dem 31/4-Theorem mit einem Wirbel (n=1) | 


A | 1 3 6 10 
ED 1 
ay 0,910 0,658 0,444 0,307 
ay 0,003 0,018 0,031 0,035 | 
a; 0,000 0,000 0,002 0,004 | 
1 dt 0,227 0,494 0,666 0,766 
20% d&g ; 
dem =s)TA <1) =e —1/4 
dicy 


o Ellipsentiige! nach ~1-Theorem (3 Aufpkt) ——_| | 
+ Ellipsentige/ nach Krienes [6] 

x Rechteckhigel nach Wiegharat [2] 
» hechteckhigel nach Blenk [1] 


oO 7 2 3 y 5 6 7 8 9 70 a 


Abb. 8. Auftriebsanstieg der ENipsen- und Rechteckflache; 
(1) Ellipsenflache als tragende Linie, GJ. (9, 4), 
(2) Ellipsenflache nach 31/4-Theorem, Gl. (9, 3), 
(3) Ellipsenflache nach 31/4-Theorem, Gl. (9, 2), 
(4) Rechteckflache nach Streifenmethode mit drei Streifen. 
(9 Aufpunkte auf einer Flachenhilfte). 


Zahlentafel 3. Koeffizienten der Zirkulationsverteilung der Streifenwirbel 


a 


| 
CI 


0 2 4 6 8 A. 


Abb. 9. Neutralpunktlage der ebenen Rechteckplatte | 
bei endlichem Seitenverhiltnis; 
(a) Streifenmethode mit 3 Streifen (9 Aufpunkte), 
(b) nach Wieghardt, 
A+: xy/l= —1/4. 


nach (5,5), Auftriebsanstieg und Neutralpunktlage der ebenen Rechteckplatte 
bei Aufteilung in drei Streifenwirbel (n=3). 


A i 3 6 10 
o 0,7151 0,4458 0,2877 0,1954 
Ti Stecifer | oe 0,0150 0,0280 0.0290 0.0264. 
. 0,0002 0,0030 0,0040 0,0041 
hi 0,1584. 0,1516 0,1075 0,0754 
2. Streifen bs —0,0097 —0,0046 0,0036 0,0069 
ie — 0.0003 —0.0021 —0.0010 0,0003 
e 0,0535 0,0685 0,0514 0,0369 
3. Streifen | cy — 0.0028 — 0.0036 0,0002 0.0024 
s 0,0001 — 0.0007 —0,0008 — 0.0002 
ib nade 
a 0,231 0,499 
sede 0,670 0,769 
at —0,179 0,228 — 0,240 —0,245 
Ca 


es | 


a te 


4 


} 


é 
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b) Rechteckflache. Die Rechteckflache als ebene Platte wurde mit drei Fouriergliedern so- 
_ wohl mit einem Wirbel (n=1) als auch mit drei Wirbeln (n=3) nach der Streifenmethode unter- 
_ sucht. Als Lésungen ergaben sich die in Zahlentafel 2 und 3 zusammengestellten Werte. Der 


Unterschied im Auftriebsanstieg und auch in der Auftriebsverteilung langs Spannweite ist hierbei 
gering. Wesentlich ist die bei drei Streifen ermittelte Abhangigkeit der Neutralpunktlage vom 
Seitenverhaltnis, die auf eine veranderte Tiefenverteilung der Zirkulation hinweist: Eine gra- e" 


 phische Auftragung gibt Abb. 9 im Vergleich mit den von K. Wieghardt ermittelten Werten, die 


hier trotz der von Wieghardt angenommenen elliptischen Zirkulationsverteilung noch in guter af 
Ubereinstimmung stehen. Der Verlauf des Auftriebsanstieges tiber A ist ebenfalls in Abb. 8 ein- ae 


a ort. 
Ca 


O27 


02 


Abb. 10. Auftriebsverteilung der ebenen Rechteckplatte iiber Spannweite. 
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10 


getragen und steht in guter Ubereinstimmung mit . 
dem von H. Blenk fiir A=6 ermittelten Wert, oa 


wahrend die Werte nach K. Wieghardt fiir groBe on 
Seitenverhaltnisse wegen der angenommenen & 
a neeh : : : : hea 
elliptischen Zirkulationsverteilung stark abwei- Ne 
chen. Auftriebsverteilung langs Spannweite und tae 
ortliche Neutralpunktlage, berechnet mit drei i 
. . NA 
Streifen, sind aus den Abb. 10 und 11 zu ersehen. x 
Ms 

. 

a 

+ by 

ang 


2% = 
i a f 


eee 


aa 


35 Ute: 
= 
Ee Fete 


as ‘ 


ae 
me 
On Oh 02) OR 0284 05% 06 07 Oo ERALO a 
Abb. 11. Ortliche Neutralpunktlage der ebenen Rechteckplatte. e 
Cat 


In gleicher Weise wie die ebene Platte wurden das Kreisbogenprofil mit «=0 sowie das nae 
S-Schlagprofil mit «,=6/4 untersucht (fir A—> 0 Birnbaumsche Normalfalle, Abb. 1). Als geo- ae 
metrische Bestimmungsgré8e des Profils sind hierbei statt des Anstellwinkels ~, die Vorderkanten- | ee 
winkel § bzw. 6 einzufithren. Die Ergebnisse der Rechnung sind in Zahlentafel 4 zusammengestellt. ee 

os 
Zahlentafel 4. Koeffizienten der Zirkulationsverteilung der Streifenwirbel nach (5, 5) i 
fiir Kreisbogen- und S-Schlagprofil bei Aufteilung in drei Streifenwirbel (n=3). sf i 
Kreisbogenprofil %g=0 S-Schlagprofil ug =0/4 tt ES 
A 1 3 6 10 1 3 6 10 3 
a, | —0,1736 | —0,0121 0,0193 0,0190 0,2583 0,1292 0,0737 0,0468 (h 
2 a | —0,0149 | —0,0192 | —0,0108 | —0,0048 0,0054 0,0129 0,0118 0,0092 i 
Streifen a; | —0,0002 | —0,0030 | —0,0031 —0,0020 Spars —0,0001 0,0011 0,0020 0,0019 . BS 
b 0,3731 0,1864 0,1093 0,0709 0,1229 0,0531 0,0296 0,0187 i 
a 6; 0,0144 0,0202 0,0165 0,0128 0,0105 0,0089 0,0056 0,0038 : 
Streifen| 5°} 9.0002 | 0,0027 | 0,0030 | _0,0025 0,0004 | 0,0019 | 0,003 | 0,009 oy 
1882 0,1060 0,0671 —0,4575 | —0,1793 | —0,0919 | —0,0553 
=: B ihe een 0,0184 0,0137 —0,0275 | —0,0334 | —0,0226 | —0,0147 
Streifen| ¢"| 9.0004 | 0,0026 | 0,0034 | 0.0029 —0,0005 | —0,0053 | —0,0052 | —0,0037 
ee Je am} o721 | 0,812 | 0,806 | 0,794 
CONE 0,300 0,543 0,698 0,785 are > ees , 
2 |= 0,733) | 0,601 |'—0,545 | —0,519 2 — - ae lids 
Ca 


er 
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Abb, 12, Hohenschichtlinien der Lastverteilung an der ebenen Rechteckplatte. 


Ferner wurde eine stetige Zirkulationsverteilung nach (7,1) ermittelt. Hieraus lassen sich 
Linien gleicher Druckdifferenz zwischen Ober- und Unterseite der Tragflache (Niveaulinien) im 
FligelgrundriB zeichnen, wie dies in Abb. 12 und 13 fir die ebene Platte und das Kreishogen- 
profil geschehen ist. Beim Kreisbogenprofil riickt der Druckberg mit abnehmendem Seiten- 
verhaltnis wegen der auftretenden negativen Saugspitze an der Vorderkante immer mehr nach 
hinten. Die Auftriebsverteilung tiber die Tiefe fir ebene Platte, Kreisbogenprofil und S- Schlag- 
profil zeigt auch Abb. 14. 

Es ist nicht schwer, durch Uberlagerung von phones Platte und Kreisbogenprofil die Werte. 
fiir ein Kreisbogenprofil mit Anstellwinkel zu erhalten. Hierbei interessiert vor allem die Null- 
auftriebsrichtung und Nullauftriebsmoment. Fir Nullauftrieb muB gelten 


sey dcq 
dug. (%g)e.-0 + ZR gh = 95 
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negative Saugspitze 


negative Saugspitze 


negative Saugspitze = 
== ee P=0 
05 
ez 
| y= 
ba 35 fas lee A=6 
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Abb. 13. Hohenschichtlinien der Lastverteilung an der Rechteckfliche mit Kreisbogenprofil. 


voraus sofort die Nullauftriebsrichtung des Kreisbogenprofils folgt. Ihre Grenzwerte lassen 


ich angeben zu B 


AOE TA Og) eager A—> 002) (0). 29 = — B- 


Das Moment bei Nullauftrieb wird dann 


Analog zum Kreisbogenprofil 1aBt sich auch jede andere Profilkontur bei sehnenparalleler An- 
tromung nach der Streifenmethode berechnen und durch Uberlagerung mit der ebenen Platte 
ler Anstellwinkel erzeugen. Auftrieb und Moment der Flache sind dann gegeben durch 


[ %_— (as) ca=0 | Cm = mg ae: Var 


we 
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Abb. 14. Auftriebsverteilung der Rechteckflache iiber die Tiefe fiir ebene Platte, Kreisbogen- und S-Schlagprofil. 


@ 


Abb. 15 zeigt das ¢my des Kreisbogenprofiles : 
in Abhangigkeit vom Seitenverhaltnis. Ebenso 
ist in*Abb. 16 die Nullauftriebsrichtung und 
die Richtung des stoBfreien Eintritts fiir das 
Kreisbogenprofil wiedergegeben. 


| 
Richtung des stobtreien Eintritts 
eal 


oe 


Nullautriebsrichtun 


2 ¥ 6 8 10 


Abb. 15, Nullauftriebsmoment der Rechteckflache Abb. 16. Anstellwinkel der Rechteckflache mit Kreicbogen pray 
mit Kreisbogenprofil; fiir Nullauftrieb und stoBfreien Eintritt; 
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Scholz: Beitrage zur Theorie der tragenden Flache. 


Zahlentafel 6. Rechenbeispiel fiir Rechteckplatte 1/b=1 (Abb. 17). 
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Tafel 6a. Gleichungssystem fiir eine ebene Rechteckplatte b/l=1 mit zwei Streifen (n=2) 
und drei Fouriergliedern (u=1, 3,5) nach. (5,8). 
a a 


Abwindgl. Abwind von Wirbel J’) Abwind von Wirbel 72 Ortlicher 
: : A 
in Aufpkt. 1. Four. Gl. | 2. Four. Gl. | 3. Four. Gl. | 1. Four. Gl. | 2. Four. Gl. | 3. Four. Gl, we 
P PISS 
(1) ay Bry > 03 3)- ay es) + by ey) + bs es) + byes} = I 
(2) ae + ase + asco + bey + bse) + obseQh = I 
(3) G83 “bagasse + ase” + bye +. 8, 63 + | byeQo = i 


PPR eusieh eds ter say, en 6. eed 10) UOT e Le" 0). “e) We: ve 


Tafel 6b. Aus Abwindtafeln zu entnehmende Abwindwerte fiir die auftretenden Koordinaten € und }. 


a rt 
o=2 o== == 
2 3 e 6 

ué ee eres SES Le ee Sey eee fe ea | keg be, 
1/2 0,1694 | —0,0533 0.018 0,1471 | —0,0121 | —0,009 0,1072 0,0191 | —0,005 
3/2 0,0774 | —0,0055 0,000 0,0733 0,0037 0,000 0,0659 | —0,0013 0,000 

€ & &3 &5 <4 3 & 2 &3 | “5 
+1/2 —0,3388 0,1066 | —0,036 | —0,2942 0,0242 0,018 | —0,2144 | —0,0382 0,010 
—1/2* 1,3388 | —3,1066 5,036 1,2942 | —0,0242 | —5,018 1,2144 6,0382 4,990 
“3/2 —0,0516 | °0,0037 0,000 | —0,0489 0,0025 0,000 | —0,0439 0,0009 0,000 
=3/2* 1,0516 | —3,0037 5,000 | 1,0489 | —0,0025 | —5,000 1,0439 5,9991 5,000 

* Die Werte der Abwindfunktionen fiir negative € ergeben sich aus der Symmetriebedingung (8, 7c). 
Tafel 6c. Koeffizienten des zu lésenden Gleichungssystems von Tafel 6a 
mit den zugehérigen Koordinaten der Aufpunkte. 
Aufpunkt von Wirbel 1°) von Wirbel J’) 

Pp eg go A en o, ao) 92 ae eo, ef, 
n/2 1,3388 | —3,1066 5,036 n/2 | —0,3388 0,1066 | —0,036 
(2) Sy r/3 1,2942 | —0,0242 | —5,018 +4] 2/3 |.—0,2942 0,0242 0,018 
(3) : 1/6 1,2144 6,0382 4,990 a/6 | —0,2144 | —0,0382 0,010 
(4) m2 1,0516 | —3,0037 5,000 m/ 2 1,3388 | —3,1066 5,036 
(5) —$ | 2/3 1,0489 | —0,0025 | —5,000 | —4| 2/3 1,2942 | —0,0242 | —5,018 
(6) /6 1,0439 5,9991 5,000 x/6 1,2144 6,0382 4,990 


Tafel 6d. Das zu lésende Gleichungssystem mit seinen Lésungen. 


Gleichungssystem: 
1,3388 a, —3,1066 a,+5,036 a, —0,3388b, +-0,10666; —0,0366;=1, 
1,2942.a, —0,0242 a, —5,018a;—0,2942b,+0,02426,+0,0186; =1, 
1,2144a,+6,0392a,+-4,990 a;—0,2144b, —0,0382b,+0,0106;=1, 
1,0516 a, —3,0037 a,+5,000 a;+1,3388b, —3,10666,+-5,0366;=1, 
1,0489 a, —0,0025 a,;—5,000 a;+1,2942b, —0,02426, —5,0186;=1, 
1,0439 a,+5,9991 a,+5,000a;+-1,2144b, + 6,0382b; +4,990b; =1. 


Aufpunkt (1): 
Aufpunkt (2): 
Aufpunkt (3): 
Aufpunkt (4): 
_Aufpunkt (5): 
Aufpunkt (6): 


Lésungen: 
a,=0,009,  a;=0,000, b,=0,124, b,=—0,007, 
T= V- b[a, sin Pa, sin3 O+a; sin 5 ] x. 
T= V- b[b, sin 9-+b, sin3 +b; sind D] x. 
Pc SES ES SE i eR 


a, =0,800, 


Institut fiir Strémungsmechanik der Technischen Hochschule Braunschweig. 


Abb. 17. Beispiel fiir die Aufstellung 
des Gleichungssystems fiir eine 
Rechteckflache A=1. 


(Eingegangen am 8. April 1949.) 
Anschrift des Verfassers: Dr. N. Scholz, (20b) Braunschweig, WodanstraBe 42. 
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Zur Ermittlung der EinfluBflachen und inneren Krifte 
umfangsgelagerter Rechteckplatten. 


Von W. Koepcke. 


1. Einleitung. Rechteckige, senkrecht zu ihrer Mittelebene belastete und an ihrém Umfang . 
laufend gestiitzte elastische Platten sind bereits in einer groBen Zahl von Arbeiten behandelt 
worden. Der tiberwiegende Teil dieser Untersuchungen erstreckte sich aber nur auf die -allseits 
drehbar gelagerte und allseits eingespannte Platte unter gleichmaBig verteilter oder dreieck- 
formiger Belastung oder einer Einzellast in der Mitte. Erst in neuerer Zeit haben Bittner! und 
Pucher? Methoden fir die Darstellung von EinfluBfeldern an rechteckigen Platten entwickelt und | 
so die Méglichkeit gegeben, fir jede beliebige Belastung die inneren Krafte zu berechnen. In 
besonderem MaBe haben die Arbeiten von Pucher dieses Problem geklart. 


Danach ist jedes EinfluBfeld fur eine beliebige statische GréBe wie z. B. Durchbiegung, — 
Biegungsmoment, Drillungsmoment oder Querkraft eine Biegeflache der Platte, die durch eine } 
gewisse, jeder dieser GréBen eigentiimliche singulare Lastgruppe erzeugt wird. Diese Singulari- | 
taten sind in dem Plattenpunkt anzubringen, fir welchen das EinfluBfeld ermittelt werden soll. | 
Sie miissen weiterhin, da sie aus Kombinationen von Einzellasten oder deren Grenziibergangen | 
bestehen, die homogene Differentialgleichung der Platte 1A w=0 erfiillen. AuBerdem haben die | 
EinfluBflachen an den Randern den Auflagerbedingungen der untersuchten Platte zu geniigen. | 


In bezug auf die Randbedingungen unterscheiden sich demnach die EinfluBflachen nicht von den 
Biegeflachen, die durch stetig verteilte Belastungen erzeugt werden. Im allgemeinen werden aber | 
weder die von den Singularitaten allein herrihrenden Biegeflachen noch die Biegeflachen fir stetig | 
verteilte Lasten, welche durch die meist leicht anzugebenden Lésungen der inhomogenen Platten- 
gleichung AA w=p(é,7)/N beschrieben werden, die Lagerbedingungen lings der vorgegebenen | 
Rander erftillen. Man schneidet deshalb aus ihnen einen die Last enthaltenden Bereich von den 
Abmessungen der zu untersuchenden Platte heraus und fiigt weitere Flachen hinzu, die unab- 
hangig von Feldbelastungen sein miissen, so dafi die Summe aller Teilflachen samtlichen Rand- 
bedingungen gentigt. Pucher hat bei Berechnung von EinfluBfeldern diese Zusatzflachen in seiner | 
ersten Veréffentlichung mit Hilfe der Differenzenrechnung, in seinem zweiten Aufsatz durch 
biharmonische Polynome gewonnen. Da aber diese Rechnungen fir jeden einzelnen Fall einen 
betrachtlichen Arbeitsaufwand erfordern, erscheint es angezeigt, ein Verfahren zu suchen, das 
eine Tabulierung zula8t und es auf diese Weise erméglicht, die standig wiederkehrenden Rechen- 
operationen auf ein Mindestmaf} einzuschranken. 


Die anschlieBenden Darlegungen stellen im wesentlichen eine Erweiterung und zahlenmafige 
Auswertung der vom Verfasser friher® ausfiihrlich entwickelten Gedankengange dar. Aus diesem 
Grunde wurden hier einzelne Abschnitte kurz gefaBt. 


2. Erlaiuterung des Verfahrens. Es ist fiir die folgenden Betrachtungen ohne Belang, ob die 
zu untersuchende Platte durch stetig verteilte Belastungen, Einzellasten oder EinfluBfeld- 
Singularitaten beansprucht wird. 

Die Gleichung der Biegeflache einer Platte mit beliebiger Belastung mu bekanntlich die 
Differentialgleichung 


_ dw 04 w 04 w p(é,n) 
AAw= yet We aera yee W (1) 
erfiillen, wobei w die Durchbiegung, p(é.7) die Lastfunktion, N die Plattensteifigkeit und &,7 
dimensionslose Koordinaten bedeuten. Bei singuléren Lasten ist p(&,7)=0 zu setzen. Die der 
Gleichung (1) geniigenden Biegeflachen, welche allein durch die EinfluBfeld-Singularitaten oder 
die partikularen Integrale der Plattengleichung bei stetig verteilten Lasten gegeben sind, haben 


* E. Bittner, Momententafeln und EinfluBflachen fiir kreuzweise bewehrte Eisenbetonplatten. Wien 1938. 
2 A. Pucher, Die MomenteneinfluBfelder rechteckiger Platten. Heft 90 der Reihe des Deutschen Ausschusses 
fiir Stahlbeton, Berlin 1938; A. Pucher, Ing. Arch. 12 (1941), S. 76. 


= 3 W. Koepcke, Umfangsgelagerte Rechteckplatten mit drehbaren und eingespannten Randern, Leipzig- 
orna 194], 


<VIII. Band 1950. Koepcke: EinfluBflachen und innere Krafte umfangsgelagerter Rechteckplatten 107 
a ON eee eee 


unachst keine Begrenzung. Man schneidet daher aus ihnen einen Bereich heraus, der die Be- 
lastung einschlieBt, und dessen Projektion gleich den Abmessungen 2a-26 der zu untersuchenden 
Platte ist (Abb. 1). Dieses so gewonnene Flachenstiick sei mit Wp bezeichnet. Es hat an allen vier 
mocdem Durchbiegungen, Verdrehungen, Kriimmungen und Querkrafte. Da im folgenden nur 
-umfangsgelagerte Platten behandelt werden sollen, wird den Querkraften keine weitere Beachtung 
geschenkt. Von den dann qn einem Rande iibrig bleibenden drei SchnittgréBen ist je nach der Art 
der Lagerung bei einem eingespannten Rand die Durchbiegung und die Verdrehung oder bei 
einem frei drehbaren Rand die Durchbiegung und Kriimmung zu beseitigen. Dies erfolgt in 
-einzelnen Schritten. Zunachst werden die Randwerte an den Ecken von wp durch den biharmoni- 
-schen Ansatz 


Oy 


é. We = a+ BE +yn + 0&y + e&? + ly? + 78 + PEN? + 168 + uni + dey + wen? (2) 


tum Verschwinden gebracht. Da an jeder der vier Ecken die Durchbiegung und die Verdrehung 

oder Kriimmung nach beiden Koordinatenrichtungen zu Null gemacht werden missen, enthalt 

der Ausdruck fiir we zwélf Freiwerte « bis i, welche so zu bestimmen sind, daB an den Ecken die 

-neugebildete Flache w)=wp+w, die endgiiltige Lagerung besitzt. Die Abb. 2 gibt ein ungefahres 

Bild von wy. Diese Flache hat am Rand r, mit Ausschlu8 der Ecken die beiden Randwerte 

Roy. To; (z- B. Durchbiegungen und Verdrehungen), die zu beseitigen sind; am Rand ry sind dies 
og? Log usw. 


Es wird nunmehr eine Flache w, gebildet, die eine Lésung der Gleichung 4A w=0 ist und des 
weiteren folgenden Bedingungen zu geniigen hat: Es miissen an den r, benachbarten Randern 
r. und r, sowie am gegentiberliegenden Rand r, die Durchbiegungen und je nach der endgiiltigen 
Auflagerung die Krimmungen oder Verdrehungen verschwinden und auBerdem in dem Ausdruck 


: 
Abb. 1. Biegeflache Wp einer rechteckigen Platte. Abb. 2. Biegeflache wy = Uy +, einer rechteckigen Platte. 
fiir w, zwei Konstante enthalten sein, die so bestimmt werden, da} 


: Ro + Ry =9, To + Ty = 0 (3) 


sind, d.h. also, daB die Randwerte R,, und T,, von w, zusammen mit den entsprechenden GréBen 
Ro, und T,, von wy, den endgiiltigen Auflagerbedingungen der gesuchten Platte langs r, gentigen. 
Die Flache w)+w, hat dann nur noch an den Randern ry, rz und ry vom Endzustand der Platte 
_abweichende Randwerte, die aber allein von w, herrihren, da entsprechend den oben aufgestellten 
Forderungen an w, 


Rye = Typ = Ryg= Tyg = Ry = Ty =9 
sind. Nun kann eine zweite Flache w, angegeben werden, die mittels zweier neuer Freiwerte die 
- Randgleichungen 
Ro tRy=0, TM+T 2 = 9 

-langs r, erfiillt und im tbrigen den durch wy+1, langs r,, rz und ry gewonnenen Zustand nicht 
-verandert, also den Beziehungen 

Roy = Tp = Ry3 = T23 = Rog = Ty, =9 
-genigt. In ganz entsprechender Weise verfahrt man mit den Randern r, und r, und erhalt endlich 
‘als Gleichung der gesuchten Biegeflache, die allen Randbedingungen geniigt, den Ausdruck 


W = Wy tw, +W2+W3 +4. 


Durch Zerlegung von unsymmetrischen Belastungen in Symmetrie und Antimetrie kann die 
Rechenarbeit zur Ermittlung von w, bis w, betrachtlich vereinfacht werden. 


Puen orden ist. Die Finder: Ton mind r, seien im Endzustand der Platte drehbar oorceal | 
dagegen kénnen r, und r, drehbar oder eingespannt sein, so da sich im ganzen die vier in Abb. : ¥ 
dargestellten Méglichkeiten ergeben. Die ee der Flache w, erfolgt nach bekannten)- 
Verfahren! und ist vom) 
a) Verfasser an anderer Stelle: i 


den. Eine Wiederholung iW 


, tre/ drehbar 


eingespannt 


a) - b) . c) af? 
Wl " é 
3 
a 
2a 
ke Tyg Tp Ty) 72 Uitig thy 
13 73 73 


Abb. 3. Rechteckplatten, bei denen beide Nachbafrander des verbogenen Randes metrie zur Ns -Achse ant f 
der Flache w, drehbar pelseert sind. gegeben werden. 
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a) Symmetrie: «) Drehbare Auflagerung von r;: 


Cof 2p — 
WwW, = yale (4.0128 }. —Besin2p 7)| cos fb € Gof hy + 
k=1,3,5.. 20% 
Gin 2p — 
ne Y, alt + En (AuGor28 5 —Becin2p-F) Joon Bs cinBy p= 2; 
k=1,3,5... eer 


6) Einspannung von r3: 


ws = A} apes + Gin2BF (AnGoi2p + = Be ein2p-£) [loos 601 ih 


| 
k= 13,9. | 


, Bul + Tap [2B s + Gof 28 (Ax Goi 28 © s Bu Gin2p-5) lleosp$ inf. 


=I, (33,05 


b) Antimetrie: ~) Drehbare Auflagerung von ry: 


of 28 
wt = eee (Af oj 28 —-— Br Gin2B §) [sim BE 180 + 
k=2,4,6 kb 
Gin 2p — 5 A 
a B14 At Go28 5 — BE Sin2p-£) sin BE cinp 3 
ae 2B; bo ( p f : 


6) Einspannung von ry 


wi = y At tf tee bated + Gin2 B= + (42 coi2p ¢ —Bt nap *)]}sinp scarey 


k=2, 4,6... 


Ay yBr rel! BE B | + Go 2p (At Gof 2h 5° — Bi ein2p- 7) || sinpéeingn. | 
k=2,4,6.. By . 4 


Die Konstanten A; und B, werden mit Hilfe einer Fourierschen Analyse der fiir die Aufloxerawal 
mafigebenden Randwerte von wo langs 7=0 berechnet. Ist wy nicht nur zur Achse 7, sondern 
auch zur Achse & (Abb. 3) symmetrisch oder antimetrisch, so kann man sich diese Tatsache 
zu Nutze machen und dadurch die obigen Gleichungen noch weiter vereinfachen. Da sich die | 
Rechnungen verhaltnismaBig einfach gestalten, erscheint hier eine Tabulierung nicht notwendig. 


' Estanave, Contribution a l’étude de l’équilibre élastique d’une plaque usw., Dissertation Paris 1900. 
* W. Koepcke, a. a. O. 
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4. Behandlung des Plattenrandes mit zwei eingespannten oder einem eingespannten und einem 
_ drehbar gelagerten Nachbarrand. Es wird wieder vorausgesetzt, daB Wy bereits bekannt ist. [m 
_ Gegensatz zu den Darlegungen im dritten Abschnitt wird w, zunachst so berechnet, als ob der 
é Rand rz unendlich weit entfernt liegt. Hat die zu untersuchende Platte in der Richtung der 7- 
| Achse die Lange 2x, so fihrt man durch w, im Abstande 4 =—2¢ einen Schnitt und entfernt auf 
\3 diese Weise den bis ins Unendliche reichenden Teil der Flache (Abb. 4). Da w,, wie im folgenden 
 gezeigt wird, langs der 7-Achse nach einer gedampften Schwingung verlauft, die rasch abklingt, 
sind die durch den Schnitt frei werdenden Randwerte bei 7=—2« klein und kénnen umso mehr 
_unbeachtet bleiben, je gréBer 2% gegeniiber 2 ist. Ist 2x wenig gréfer, gleich oder kleiner als 2, 
| "so wird es erforderlich sein, eine zweite Flache w; auf entsprechende Weise wie w, zu bestimmen, 
die die Aufgabe hat, die 
| kleinen Randwerte von 
4 w, langs n=—2« zu be- 
) seitigen. Auch wi ist ein 
-unendlich langer Halb- 
-streifen, der am Rande 
_r, durchzuschneiden ist. 
| Sind die durch diesen 
| Schnitt — entstandenen 
oe Randwerte noch so groB, 
-daB sie sich bemerkbar 
- machen, so ist eine dritte 
Flache wi’ zu berechnen, 
| was aber nur in ganz sel- 


fref drehbar 


elngespannt 


; z : ™% iKige Qe) Ty\ Te OPS Ui "Y 
tenen Fallen notwendig 20 
-sein wird. Die Einzel- 
heiten, die bei Bestim- pers 
1 


H 3 HT i 73 
H Ht i } 


iV 
| t 
| ! i " 


mung der Flachen w; bzw. 


ZA . 
peu beachten sind, lak Abb. 4, Rechteckplatten, bei denen beide oder ein Nachbarrand des verbogenen Randes 
den weiter unten ausfihr- der Fliche w, eingespannt ist. 


licher besprochen. 


Durch die Benutzung eines Halbstreifens geht das Seitenverhaltnis der Platte in die Rech- 
nungen nicht unmittelbar ein, und man gewinnt daher die Moéglichkeit einer weitgehenden 


Tabulierung. 


Im folgenden werden die theoretischen Entwicklungen fiir den Fall dieses Abschnittes dar- 
' gelegt, w, bei einem oder zwei eingespannten Nachbarrandern von r, zu berechnen. 


Es wird vorausgesetzt, da die Randwerte von wy langs r; in Symmetrie- und Antimetrie- 
| teil, bezogen auf die Achse 7, zerlegt worden sind. 


a) Symmetrieteil. Fur w{ wird der Ansatz} 


Ww; = » Cye?sk” (@sk& sin sh & + sin? wsk COS Wsk&) + | 

E 

rs vk & (4) 
+ » Cy?" (sk Esin Wsk& + sin? Wsk COS sk &) | 

kK 
‘® gemacht. Bei 
. On >= Pru = -- 2,106196 ys Dl 25300%0% 
sy = Po + Got = + 5,356269 + 1,551575 1, 
Os3 = P3 + 93t = + 8,536683 + 1,7755447 usw. 


, verschwinden langs = + | die Durchbiegungen und Verdrehungen von w;. Die mit einem Quer- 
| strich versehenen GréBen sind konjugiert komplex zu denen ohne Querstrich. Es bedeuten 


Ce = ReC, LiQmC, und Cy = ReC,—i9mC,. Fuhrt man noch 


®;(E,) = Re Ox (E, n) +1 9m Oi(E, 7) = Ag 7 (ask € sin Wsk € + sin” Wsk C08 Wek &), 
Di (E,n) = Re Di (E,n) —i 3m Dx (E, y) = eS" (Gen € sin Ost § + sin? Wsk COS sk E) 


1 W. Koepcke, a. a. O. 
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ein, so erhalt man ; 


wi? = D)[Ce@x(é,) + Ce Pelé, m)] - 
B 
Nach Ausmultiplikation ergibt sich 


ws — 2 >' [me Ci Re De (E,) + Sm Cy Jm De (E, 79) | , = (5) ye 
: | 


und man erkennt, daB durch die Einfihrung der konjugiert komplexen Groen der Ausdruck fur | 
w’ reell geworden ist. Die Gleichung (5) enthalt die beiden voneinander unabhangigen Konstanten | q 
Re C, und Sm Cy, die bendtigt werden, um die Randbedingungen langs r, zu erfiillen. | 


a) Es werde zunachst angenommen, da r, im Endzustand eingespannt ist, wi also bestimmte 
Durchbiegungen und Verdrehungen langs 7 = 0 anzunehmen hat, und zwar so, dab 


wit f=0, (>) +y.=0 (6) | 


werden, wenn man unter wi, und (dw}/d7), die Durchbiegungen und Verdrehungen von w} und |) 
unter f, und gs die Durchbiegungen und Verdrehungen von wy am Rande r, versteht. Setzt man 
zur Abkiirzung 


0 ws 


0% 


PD, = Re Dy +i Sm Ox = asx E sinasz & + sin? Wsk COS Wsk €, 
Dy, = Re Oy —i 3m Oy = Wek € sin Wex E + sin? Wsk COS Wsk € , 


so werden bei r, fiir 7 = 0 


Wi » (Cy ®, + Cr Dx)» 
rs 


doy : im | ay 
( j ) = >) (Chas De + Cesk Dr) « 

1] 1 k | 
Es ist far die Ubersichtlichkeit der Rechnung zweckmaBiger, im folgenden im Komplexen | 


weiterzuarbeiten. Die noch unbekannten Konstanten C, und Cx werden aus einer Minimal- | ‘ 
bedingung berechnet, und zwar wird die Forderung aufgestellt, da die Summe der Quadrate der 
Abweichungen der zu berechnenden Randwerte von den gegebenen ein Minimum wird. Es ergibt 
sich dann 
+1 : 
3 = fi iLf. at >) (Ci. Oy + Cy D,) | —E [ys + > (Cr ose Dy + Cesk d)) | d& = Min. (8) 
4 k ke 

Dieser Minimalansatz ist in der angeschriebenen Form ohne weiteres statthaft, da der Ausdruck (8) 
reell ist. ; 

Man bezeichnet nunmehr mit k die laufende, mit h eine bestimmte Ordnungszahl der Reihen- 
glieder in (8) und bildet die Ableitungen nach den beiden Unabhangigen Ste C, und 3m Ch. 


Oe 0G Ea a Ge eS a tar RO 20 ul a0, 
aC, ame Cy aC, agm, ” 70, oReG, 9G, OMG,” 


Setzt man die so gewonnenen Ausdricke gleich Null, so erhalt man 


+1 : 
ie ie Shr 2 (C.D 4+ Ci, Dr) | D, 92 [ps ie ye (Cy sk PD; aL Cr@sk Dr) | Osh Pr} d& = 0 ’ ; 
nan, (9) | 

| 2 ne Si. ne (CDi. + CrDx) | On + 2 [ys ++ S (Chosk Di + Cr Ose Dx) | wn Dy dé = 0. 
= 


= 
k ke J 


Die Trennung der Gleichungen (9) in Real- und Imaginarteil ergibt ein Gleichungssystem fir die 
gesuchten Konstanten ReC, und SmCp. Es ist aber nicht immer angebracht, diesen Weg zur 
Ermittlung von Xe C, und $mC; einzuschlagen; denn der Minimalansatz soll erreichen, da das 
Integral tiber jede einzelne der eckigen Klammern in den Gleichungen (9) in der Grenze zu Null 
wird. 

Die Ansatze (9) kénnen aber auch verschwinden, wenn z.B. das erste, durch f, gekennzeichnete 
Glied eines Integrals einen groBen Wert hat und das zweite, durch gs gekennzeichnete ebenso 
grof} ist wie das erste Glied, aber ein umgekehrtes Vorzeichen besitzt. Unter solchen Umstanden 
ware die Minimalforderung nicht erfiillt, und es wiirden schlechte Ergebnisse zu erwarten sein. 


XVILI. Band 1950. Koepceke: EinfluBflachen und innere Krifte umfangsgelagerter Rechteckplatten. ~ i11 
LLL 


Man kommt dagegen stets zum Ziel, wenn man beide Integrale (9) addiert. Es ergibt sich dann: 
E ea. 

4 hee. (Cx Dy + CPx) | Re By + 

; 2 ee i 

ea + [e+ a (Cx sk Px + Cr@sk Px) | (prRe On— gus Pr) | dic x0 

Fa k 


_ Der Minimalansatz hat nach dem oben Gesagten nur einen Sinn, wenn in der letzten Gleichung 
jede der eckigen Klammern wiederum fiir sich Null wird, also wenn 


2 +1 
s [LG+D (GO. + Cy] Reade = 0, (10a) 
; k 
= 
+1 
[le + >) (Cheose De + Cu@st Br) | (pa Re On— qu 3mDy) dE = 0 (10b) 


% = 
ist. 
Die Gleichung (10b) kann aber auch unbrauchbar werden, namlich wenn pa Re Op —qn Sm OD, = 0 
ist. Um diese Méglichkeit auszuschlieBen, multipliziert man die runde Klammer aus und erhalt 
+1 
iD [ys Se oy: (Cysk Dy + Cr Wsk P;) | pave DB, dé — 
Te 
ai 
+1 ty) 
— f [os + >» (Crosk Px + Cesk Pr) | qu SmOnrdé = 0. 
ke 


—1 


_ Auch hier ist der Minimalansatz nur erfullt, wenn jedes einzelne dieser beiden Integrale fir sich 
| Null ist. Da auBer (10a) noch eine zweite Gleichung zur Ermittlung der Konstanten SteC, und 
| SmC, bendtigt wird, verwendet man weiter nach Division durch py» das erste Integral in (11) und 
| erhalt als endgiltige Bestimmungsgleichungen 

pl 


[Us a » (Cx Op + Cr Dx) | Re OpxdE = 0, 
E 


+1 
Ft [ps = 3 (Crcosk Pr + Ch Osx Pr) | Me Did E220. 
k 
Si 


Schreibt man ftir Cz, Crs Dr, DP; die Real- und Imaginarteile, so entstehen nach kurzer Zwi- 
) schenrechnung die Gleichungen 


tei Ber 41 
i 
> [ste Ce [ xe Dy. Re Ord E + 3m Ce [8m Op Ke Prd E] 25 fs Re Ord E = Zip, 
id ay ‘ 1 hq 
Ste ary 
>) |e CL [ Se Ddte Prd E— gp [sm Ke Orde] + (12) 
Ss = 
ae - 41 ; apt! 
+ Sm Cr [ ge [Re Pre Pid é “+ pk [3m Or ve ide] eee fp Re Duds = Zing. 
asi —l —1 


W Auf den linken Seiten durchlaufen k und h die Zahlen 1, 2,3...usw..k Sh. Fur die Genauigkeit 
8) der Rechnung ist es in den meisten Fallen ausreichend, bis h = 2:zu gehen, also vier Gleichungen 
fiir die vier Unbekannten ReC,, 3mC,, Ke Cy, SmC, zu lésen. Bei hoheren Anspriichen kann auch 
} noch h=3 hinzugenommen werden. 


| B) Ist der Rand r, im Endzustand nicht eingespannt, sondern drehbar gelagert, so hat wi 
ff langs. 7 =0 vorgegebene Durchhbiegungen (—fs) und Kriimmungen (—ys) anzunehmen, also 
Q* 
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entsprechend (6) die Gleichungen 


4 


0? ws 


tee Naan ame 


zu erfillen. Hierin bedeutet ys die Kriimmung der Flache wy, langs des Randes r,. Wahrend sich 
ws, gegentber (7) naturgemaB nicht andert, wird 


( e : a »: (Cy, 0%, D, + Ci, OD) r 


0° Kk 


Ein Vergleich dieses Ausdruckes mit (7) zeigt, daB an Stelle von wsx und @s, ihre Quadrate | 


treten, im tbrigen aber keine Veranderung vor sich geht. Setzt man wy), =s,+t,t bzw. 


>, = s,—t,i, so kénnen sdmtliche Ableitungen des Abschnittes ~) tbernommen werden. In 


+1 


+1 
(12) sind nur p, q, @s und Zig=—1 f gyskeD, dé durch s,t,y. und Ziy=—2. {pote @, dé zu ersetzen. — 
—1 —1 | 


b) Adtimetrieteil, Ganz ahnlich wie bei dem Symmetrieteil gestalten sich auch die Rech- | 
nungen fir den Fall der Antimetrie. Dabei ist 


wt = >) Ch ek” (aa Ecos War — cos® warsin Waré) + 
k 
2 (13) 
=k > Ck e°2k” (Wan E COS Mak & — COS? Wak Sin Wak €). 
k 


Die GréBen wax sind hier 
. Oa =pi tqri=+ 3,748838 + 1,3843397, 
Wag = p2 + q2t =+ 6,949980 + 1,6761051, 
Wag = ps + q3t = £10,119259 4+ 1,8583847. 


+ 1 die Durchbiegungen und Verdrehungen | 


Fur diese Werte verschwinden an den Randern & = 4 


der Flache (13). i 
Setzt man 


We(En) = Re Va (E,7) +i 3m Ve (Ey) = ek" (War € cos War € — Cos? Wak Sin Wak €) , 


(En) = Re Ve (E,) —i 3m Va(E, 7) = eek” (Dak € cos Dak & — cos? Wak sin Dak €), 


so ergibt sich entsprechend den Entwicklungen beim Symmetrieteil 
wt = 2 pa [Cz Pu(é,) + Ct (En) ] 
und nach Ausmultiplikation 
w? =2 >, [Re Cf Re Pe(E.m) + SmCf Sm Vi (E, n)|- (14) 
Die Ermittlung der Konstanten ReCf und SmCf erfolgt ebenfalls wie beim Symmetrieteil. — 


a) Ist der Rand r, im Endzustand eingespannt, so erhalt man analog zu (12) 


S511 +1 +1 
> [me ce f xeWere Pade + gmcz f gm Pade Pr de] eee J foe Wa dé = Zky, 
k 
Ni) ata €. 
+1 +1 
D) {rect [pif te We te Pads —gt f gm Pave Ys de] + (15) | 
=i =i 


+1 +1 


oe] 
+ gmt [gt [ ste Paste Ps dé +pt [ gm ¥,Re Ps de]| == f[ vote Pi dé = Zh, 
a ai 24 


B) Im Falle der drehbaren Auflagerung des Randes r, sind entsprechend den Darlegungen fir 
+1 41 
den Symmetrieteil p*, q*, ga, Zig= — aah ga ke Yi, dE durch s*, t*, pa Liy = Sy Re Wn d& zu 
=k —-1l 


ersetzen, wobei s* und t* die Real- und Imaginareteile von w; darstellen. 
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L Aus der Antimetrieflache w{ fiir Einspannung beider Nachbarrander ergibt sich auch die 
_ Flache w, bei drehbarer Auflagerung des einen und EKinspannung des anderen Nachbarrandes. 
_ Die Schnitte 7=konst durch w? haben im wesentlichen den in 
_ Abb. 5 dargestellten Verlauf. Betrachtet man nur den rechten, 
_ ausgezogenen Teil der Kurve, so erkennt man, daB an der Stelle 
_ €=0 Durchbiegung und Kriimmung verschwinden, wogegen bei 
_ €=1 Durchbiegung und Verdrehung zu Null werden. Demnach 
_ erfillt also die Antimetrieflache w? bei €=0 und €=1 die Be- 
_ dingungen der Abb. 4e bis 4h an den Nachbarrandern von rs 
_ und man kann alle Ergebnisse des Abschnittes 4b) fir 0 S&< 1 
 tbernehmen. 


Abb. 5. Schnitt 7 =konst durch eine 
antimetrische Flache Ce 


c) Anstatt die Gleichungen (9) zu addieren, kénnte man sie auch voneinander subtrahieren. 
_ Die Entwicklungen wirden dann ganz entsprechend wie oben verlaufen und fir Symmetrie 
_ und Antimetrie Gleichungssysteme liefern, die wie (12) bzw. (15) ausséhen, nur waren Re@®, durch 
 §m@rz bzw. Re Y, durch $m Y, zu ersetzen. Es verlohnt sich nicht, neben (12) und (15) auch noch 
diese neuen Gleichungen numerisch auszuwerten, da Zahlenrechnungen zeigen, daB bei ihrer 
_ Verwendung keine Verbesserungen der Ergebnisse erreicht werden. 


5. Erlauterungen zu den Zahlentafeln. Die Verwendung eines Halbstreifens bei der Ermittlung 
- von w, schaltet, wie bereits bemerkt, das Seitenverhaltnis der zu untersuchenden Platte aus und 
erméglicht so die Aufstellung von Zahlentafeln, die den Umfang der Rechnungen auf ein geringes 
MaB vermindern. 

Zunachst kénnen die Gleichungen (12) und (15) fiir Symmetrie und Antimetrie bei vorgege- 
benen Durchbiegungen und Verdrehungen oder vorgegebenen Durchbiegungen und Kriimmungen 
am Rande r, ein fiir alle mal aufgestellt und in Abhangigkeit von den Z-Werten gelist werden. 
Die Lésungsmatrizen fiir die Konstanten Re C;, und 9m C; bzw. Re Cz; und SmC; stehen im Tabellen- 
| teil unter Ab und Bb. Es sind jedesmal drei Falle unterschieden, namlich h=1 (erste Naherung), 

 h=1, 2 (zweite Naherung) und h=1, 2, 3 (dritte Naherung). Je héher der Grad der Naherung ist, 
desto genauer sind die Ergebnisse. Bei normalen Ansprichen geniigt es im allgemeinen, sich der 
zweiten Naherung zu bedienen. 

In den Lésungsmatrizen sind nur die GréBen Zny, Zig, Ziy baw. Zigs Lig» Ziy von der jeweils 
gegebenen Flache w, abhangig und in jedem Einzelfall besonders zu berechnen. Da die geschlos- 
sene Auflésung der Integrale Z, der Gleichungen (12) und (15) meist sehr zeitraubend oder mit- 
unter gar nicht durchfihrbar ist, wendet man am besten stets die Siumpsonsche Regel an. Aus 
diesem Grunde sind im Tabellenteil unter Aa und Ba die Funktionen Re®,, ReD, baw. Re VY, 
Re Y, in Intervallen von 0,05, die Funktionen Re®,, Re Y, in Intervallen von 0,025 langs der 
&-Achse zwischen €=0 und =] tabuliert. 

Sind die Konstanten eC, und 3mC; bzw. KeCf und SmCfz bekannt, so setzt man sie in 
(5) oder (14) ein und gewinnt auf diese Weise die gesuchten Biegeflachen wi, und w{,, fir h=1 
bis h=3. Die zugehérigen Zahlentafeln sind unter Ac und Be zu finden. Die Flachen wi}, und 
wi, sind fir 0S € [1 und 7 [0 in Intervallen von 0,2 nach beiden Achsrichtungen berechnet. 
Die Intervallteilung der Flachen w{,, ergibt bei ihrer Anwendung auf die in Abb. 4e bis 4h dar- 
gestellten Platten ein Maschennetz, das doppelt so weit ist als dasjenige der Platten der Abb. 4a 
bis 4d. Aus Raummangel konnten die Tafeln fiir eine engere Teilung hier nicht zum Abdruck 
kommen, doch wird bei den meisten Aufgaben der Platten der Abb. 4e bis 4h auch die groBere 
Maschenweite ausreichen. 

Des weiteren enthalten die Zahlentafeln unter Ad und Bd Werte fiir die Verdrehungen der 
Flachen w,,) in Richtung der 7-Achse, die bendtigt werden, wenn die Verdrehungen am gegen- 
aberliegenden Rande r, anzugeben sind. SchlieBlich sind in den Tafeln Ae und Be die Krim- 
mungen von wi, und wi, nach 7 und in den Tafeln Af und Bf nach € enthalten, mit deren. Hilfe 
die Biegungsmomente der Flachen w, berechnet werden kénnen. 


Es sind nun noch einige Bemerkungen tiber die Bestimmung der Flachen w, und wy’ 
~zu machen, woriber schon oben kurz gesprochen wurde. Der unendliche Halbstreifen der 
Flache w, wird am Rande r, durchgeschnitten, es entstehen daher an dieser Stelle freie Randwerte, 
die durch geeignete MaBnahmen beseitigt werden miissen. Das geschieht durch eine Flache uj, 
die vom Rande rz beginnend nach r, zu abklingt und bei rz die umgekehrten Randwerte von w, 
hat. 
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an ae 
1 1 
‘5 Zig =— = f g-Re nde, Zry = — = f poe nde, 
aah) soll 


é KeD, Re®D, Ke D, 
a EEUU EEE EEE EEE EEEEEEEEEEEEE? 
0 1,650614 2,057848 2,278461 
0,025 2,192129 
0,05 1,639815 1,946111 1,9387.92 
0,075 1535040 
0,10 1,607704 1,622726 1,007345 
0,125 0,390423 
0,15 1,555121 1,121842 —0,275088 
0,175 —0,945069 
0,20 1,483452 0,497052 —1,574865 
0,225 —2,136298 
0,25 1,413542 —0,189140 —2,549398 
0,275 — 2,826858 
0,30 “ 1,290879 —0,863500 — 2,933954 
6,325 —2,863187 
0,35 1,175082 —1,459256 —2,609374 
0,375 — 2192353 
0,40 1,050280 —1,919314 —1,633351 
0,425 —0,964913 
0,45 0,919802 —2,203729 —0,226449 
0,475 0538366 
0,50 0,787139 —2,294459 1,284485 
0,525 1,968657 
0,55 0,655836 —2,197245 2,551824 
0,575 3,002869 
0,60 0.529391 —1,940726 . 3,299649 
0,625 3,431123 
0,65 0.411123 —1,572383 3,397975 
0,675 3,212388 
0,70 0,304072 —1,151644 2,897191 
0,725 2483859 
0,75 0,210860 —0,740983 2,010014 
0,775 1,516209 
0,80 0,133574 —0,395624 1,042351 
0,825 0624057 
0,85 0,073646 —0,153452 0,289071 
0,875 0,054271 
0,90 0,031733 —0,026188 —0,076746 
0,925 —0,115062 
0,95 0,007580 +0,006693 —0,086643 
0,975 —0,031239 
1,00 0 0 0 


Dritte Ndherung. 
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Tafel Ab. Symmetrie. 
I, Einspannung des Randes 1. 
Erste Néherung. Zweite Néherung. 
| Lif Lio Lif Lio Zof Zoo 
1 —0,269574 | +0,215644 ec, —0,341522 | +0,229421 | —0,283047 | +0,024799 
SmC, |.—0,588194 | +0,164048 Sea —0,600108 | +0,165184 | —0,040241 | —0,004027 
| 
Re C, —0,067188 | +0,018486 | —0,296860 | +0,063074 
gm C, —0,064528 | +0,013387 | —0,259818 | +0,029559 
Dritte Néherung. 
| Zif Zio | Zs Zoo Zsf Zs 
Rie C, —0,315971 | +0,226601 | —0,161336 | +0,010498 | +0,165547 | —0,012802 
ym C, —0,586045  +0,163722 | +0,026479 | —0,011054 | +0,082546 | —0,004084 
Je C; —0,094988  +0,021272 | —0,428442 +0,075994. | —0,153336 | +0,004672 
Sm C, —0,053074 +0,011976 | —0,204819 | +0,021773 | +0,088167 | —0,010563 
Ses —0,028971 | +0,003591 | —0,139166 | +0,019892 | —0,225034 | +0,027423 
Sm C, —0,024555 | +0,002682 | —0,116877 +0,013487 | —0,156497 | +0,011409 
II. Drehbare Auflagerung des Randes ry. 
Erste Ndaherung. Zweite Ndaherung. 
| ee | Ziy Zif Ziy Zof Z2y 
Re C, +0,022352 | +0,051193 He Cy —0,041319 | +0,054026 | —0,289619 | +0,000881 
oi, —0,366115 | +0,038944 Snrc, —0,391096 | +0,039746 | —0,111200 | —0,001432 
: Je C, —0,023083 | +0,001991 | —0,112546 | +0,005838 
sm C, —0,037844 | +0,002065 | —0,175123 | +0,002706 


a aa al BS A UN I el eh a 
| Lif | Ly Zof Zoey Zs | Zs y 
Kec, —0,009566 | +0,053721 —0,140673 —0,001061 +. 0,223200 —0,000767 
gmc, —0,368772 | +0,039542  —0,006364 —0,002642 +0,152440 —0,000122 
Re C, —0,043606  +0,002170 | —0,209003 | +0,006844 —0,137067 —0,000174 
astC, —0,033118 | +0,002002 | —0,153137 +0,002170 | +0,040365 —0,000779 
Jte C, —0,012214 | +0,000151 | —0,056943 | +0,001219 —0,099485 +0,001563 
Ym C; —0,014497 | +0,000146 | —0,067928 +0,000986 —0,104787 +0,000619 
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r9o— 0 0 0 0 0 0 0 0 0 0 0 0 (pea 
Ce 0 0 0 0 0 0 0 0 10000°0— 0 10000°0— 0 0 
0:9 0 0 0 0 0 0 10000°0— 0 10000°0— 0 10000°0— 0 0°9— 
oS 0 0 ‘ 0 0 0 0 T0000°0— 0 10000°0— 0 10000‘0— | 10000°0+ | 8°¢— 
oS 0 0 0 0 10000°0— 0 10000‘0— | 10000°0+ | 10000‘0— | 10000°0+ | 20000°0— | 10000‘0+ | 9°¢— 
nS 0 0 0 0 10000°0— | 10000‘0+ | 10000‘0— | 20000‘0+ | 20000‘0— | 20000°0+ | 20000°0— | 20000°0+ | r°S— 
(B35 —— 0 0 0 0 10000‘0— | 10000‘0-+ | 10000‘0— | z0000‘0+ | z0000‘0— | g0000‘0+ | Z0000°0— | €0000°0+ | Z°S— 
v= 0 0 0 10000‘0+ | 10000‘0— | ¢0000‘0+ | 10000‘0— | S0000‘0+ | 10000‘0— | L0000‘0+ | 10000‘0— | L0000‘0+ | 0*s— 
ie 0 0 0 T0000°0+ | 10000‘0— | ¥0000‘0+ 0 L0000‘0-- 0 01000°0+ 0 11000°0-+- | 8°F— 
oP 0 0 0 Z0000'0+ | 10000‘0+ | 90000‘0+ | Zo000‘0+ | T1000‘0+ | ¢0000‘0+ | S1000‘0+ | F0000°0+ | 91000‘0+ | 9°F— 
ry 0 0 10000°0+ | ¢0000°0-+ | €0000‘0+ | 60000°0+ | 90000‘0+ | 91000‘0+ | 60000‘0+ | 12000‘0+ | T1000°0+ | €2000°0+ | ¥F— 
Gila 0 0 Z0000°0+ | ¥0000‘0+ | 20000‘0+ | Z1000‘0+ | St000‘0+ | 12000‘0+ | 12000‘0+ | 8z000‘0+ | €2000°0+ | 0€000‘0+ | 3°r— 
0°r— 0 0 £0000°0+ | S0000‘0+ | S1000‘0+ | 91000°0+ | 62000‘0+ | 92000‘0+ | T¥000‘0+ | HE000‘0+ | SP000.O+ | LE000°0+ | 0°F— 
8°E— 0 0 L0000‘0+ | 90000‘0+ | 82000‘0+ | 81000°0+ | z2S000°0+ | 62000‘0+ | 22000‘0+ | L¢E000‘0+ |} 62000°0+ | 6€000°0+ | 8°¢— 
9°E— 0 0 F1000‘0+ | 90000‘0+ | 8F000‘0+ | Z1000‘0+ | 88000‘0+ | 92000°0+ | 61100‘0+ | 0¢000°0+ | TE100°0+ | 1E000°0+ |} 9°E— 
Weer 0 0 Z2000°0+ | S0000‘0+ | LL000°0+ | 01000°0+ | 6E100°0+ | 80000‘0+ | L8100‘°0+ | F0000°0+ | $0Z00°0+ | 20000°0+ | FE— 
Ge 0 0 Se000'o+ | . 0 LIL00‘0+ | Z1000‘0— | 60200‘°0+ | Sgo00‘0— | 62200°0+ | 8S000‘0— | S0g00‘o+ | 29000‘0— | z‘E— 
Once 0 0 ZS000°0+ | Z1000‘0— | 0L100°0+ | 8S000‘0— | 66200‘0+ | €2100°0— | S6E00°0+ | O8LTO0‘O— | OgFO0‘O+ | €0Z00°0— | 0°E— 
SiGe 0 0 €L000°0+ | 9€000‘0— | 7€z00°0+ | grT00‘0— | z0F00‘0+ | $8z00°0— | 9zS00‘0+ | g0F00‘0— | 02S00‘0+ | 8PP00‘O— | 8°t— 
oo= 0 0 S6000°0+ | 8L000°0— | 96200°0+ | 26200‘0— | 10S00‘0+ | 09S00°0— | ¥F900‘0+ | LLLOO‘O— | ¥6900°0+ | 6S800°0— | 9°%— 
TaGes 0 0 SI100‘0+ | 6F100°0O— | T¥E00‘0+ | PESOO‘O— | ESS00‘0+ | ZOOTO‘O— | 26900°0+ | ELELO‘O— | 6EL00°O+ | EIS1O°O— | H°2— 
Ot6= 0 0 IZ100°0+ 09Z00°0— SZE00'O+ | 01600°0— | T8P00‘0+ | 9L9T0‘O— | 19S00‘0+ | FLZZ0‘0— | F8S00°0+ | 667Z0°0— | Zt— 
Ores 0 0 16000°0+ | LZ100‘0— | FLI00°0+ | Z9F10°O— | SFIO0‘O+ | €S9Z0‘0— | 09000°0+ | 69S¢0°0— L1000°0+ | O16€0°0— | 0°¢— 
ot 0 0 01000°0— | $9900°0— | 8PrEgo0‘o— | 0€zz0‘0— | 28900‘0— | 066€0‘0— | PZTLO‘O— | ZzEso‘o— | POETO‘O— | €18S0‘°0— | 8° I— 
OF = 0 0 ZZ00°0O— | $8600°0— | FITLO‘O— | 1EzEg0‘O— | BLEZO‘O— | £69S0°0— | O8FSO‘O— | ZZSL0°0— | FI6E0‘O— | Z6180‘0— | 9°T— 
7 i 0 0 00L00°0— | O8E10‘0— | FSLZ0‘0— | 1zPPO‘O— | 89FSO‘O— |. 8F9L0‘O— | 8ZLLO‘O— | 16660°0— | 00980‘0— | LE8OT‘O— | FI— 
Gli= 0 0 00S10°O— | FI8LO‘O— } F6SS0‘O— |! TE9S0°O0— | FZLOL*O— | 00S60°0— | 9L8FL‘O— | ZIZZL‘O— | O9F9T‘O— | SLIET‘O— | 21I— 
= 0 0 S¥8Z0°0— | S8IZ0°0— | PZZOL‘O— | 29F90°O— | O9T6T‘O— | S9POT‘O— | 1SZ9Z°0— | L80ET°O— | 62687°0— | 696E1T°0— | O°I— 
Si Oi 0 0 896F0'0— | 8LZZ0°0— | S8ELI‘O— | IZ190°O0— | 610zE‘0— | 82060°0— | ZHPEv'O— | ZOSOT‘O— | SZLLP‘O— | TI60T‘0— | 8°0— 
9°0— 0 0 SST80°0— | P69T0‘O— | S06zLz‘O— | O8TEO*O— | ZE90S*0O— | 1zSz0‘0— | LTOI89‘0— | 1Z800°0— | 86ShL‘0O— | 72000°0+ | 9*°0— 
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Tafel Ba. Antimetrie. 


41 ey 41 
1 ee 1 ss 1 
Lig = —= | fate Padé, Lip = — = [ pete Padé, Zity=— = | vote Pad €. 
Ly —1 


—l 


& Rev, Ste VY, Fie Y%, 
0 0 0 0 
0,025 —0,800541 
0,05 —0,421672 —0,961728 —1,523013 
. » 0,075 —2,096286 
0,10 —0,822079 —1,753142 —2,462502 
0,125 —2,582301 
0,15 —1,181172 = 2,233304 —2,438405 
0,175 —2,037317 
0,20 —1,481925 —= 25315356 —1,408971 
0,225 —0,604275 
0,25 —1,710324 —1,980976 0,309129 
0,275 1,243558 
} 0,30 —1,857360 — 1,283930 2,143543 
a 0,325 2,902054 
i 0,35 —1,919035 —0,340019 3,458610 
| 0,375 3,759855 
a 0,40 —1,896770 0,693793 3,773599 
% 0,425 3,491808 
0,45 —1,797316 1,647925 2,931510 
i ogi 2,133536 
0,50 —1,632287 Qe lote 1,159058 
0,525 0,084140 
0,55 —1,417332 2,759864 —1,006881 
0,578 —2,029062 
sy 0,60 —1,171001 2,775293 — 2,904503 
“Si 0,625 ; —3,569649 
iy 0,65 —0,913350 2,452408 —3,981040 
0,675 : —4,119370 
0,70 —0,664202 1,890598 —3,991171 
0,725 —3,627918 
0,75 —0,441357 1,231276 —3,082498 
i 0,775 —2,423431 
“ 0,80 —0,259325 0,624315 —1,727105 
0,825 —1,068772 
0,85 —0,126455 0,189557 —0,514006 
0,875 —0,108841 
0,90 —0,044630 —0,018997 0,129783 . 
0,925 0,199451 
0,95 —0,007518 —0,037757 0,152789 
0,975 0,055540 
1,00 0 0 0 


Erste Ndherung. 


Tafel Bb. Antimetrie. 


I, Einspannung des Randes r,. 


Zweite Naherung. 


Re CT —0,313972 | +0,104799 Rect —0,426773 | +0,119591 | —0,213430 | +.0,011169 
SmCfT | —0,358745 | +0,056996 Sm Cy —0,339456 | +0,052612 | +0,051648 | —0,010545 
ste Cp: —0,111536 | +0,021033 | —0,263385 | +0,040878 
Sm CF —0,098347 | +0,014630 | —0,200239 | +0,017807 
Dritte Naherung. 
ce et, is zy me 5 ay 
| ’ 
E ect | —0,365116 | +0,114063 | —0,039227 | —0,004932 | +0,171664 | —0,009583 
gm Cy —0,322514 | +0,051280 | +0,097397 | —0,013886 | +0,034350 | +0,000239 
Re CT —0,159371 +0,024963 | —0,394478 | +0,051292 | —0,108614 | +0,001930 
Sm Cy —0,069309 | +0,011729 | —0,114835 | +0,008504 | +0,101196 | —0,009073 
Re CF —0,050345 | +0,005243 | —0,150488 | +0,017365 | —0,190072 | +0,019005 
nics —0,042914 | +0,003912 | —0,121970 | +0,011577 | —0,123865 | +0,007653 
II. Drehbare Auflagerung des Randes 1. 
Erste Naherung. Zweite Naherung. 
eee aes Zn, ie ay Bs 
Re Ct —0,090748 | +0,013977 Re CT —0,177935 | +0,015641 | —0,195511 | —0,000017 
gm CF | —0,237342 | +0,007602 Siu Ore —0,239881 | +0,007326 | —0,000059 | —0,001155 
=a Re Cy —0,042857 | +0,001629 | —0,110233 | -+0,002887 
Sin CF —0,057378 | +0,001440 | —0,134671 | +0,001219 


a 


Dritte Naherung. 


ST ne Saat IESsEEEESCLEC ARS 
| 


7 ae - he Be iy: 

Ste CT | —0,119367 +0,015220 —0,034302 —0,001354 +0,179424 —0,000310 
Sm Ct —0,216370 +0,007182 | +0,064172 —0,001561 +0,065657 +0,000162 
Re CF —0,077003 +0,001847 | —0,203672 +0,003519 —0,097443 —0,000141 
Sm Cz —0,042392 +0,001295 | —0,092700 +0,000682 +0,.055360 —000,0504 
Re CF —0,022705 +0,000231 —0,063828 +0,000880 —0,087079 -+0,000908 

: Sm C3 —0,025437 +0,000201 —0,070367 +0,000676 —0,082568 +0,000343 


a 


Ingenieur-Archiy the 
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Besonders einfach gestaltet sich die Ermittlung von w;, wenn r; drehbar aufgelagert werden y 
soll. Die betrachtete Platte habe das Seitenverhaltnis a/b = 2/2. Bei n=—2a hat w, bestimmte | 


Durchbiegungen w,, und Kriimmungen (0720,/07?)3,.siehe Abb. 6. Der durch den Schnitt bei rs 
abgeléste Teil von w, wird um 180° um eine Achse gedreht, die senkrecht zur (€,7)-Ebene im 


Punkt €=0, 7 =—2« steht und, wie in Abb. 6 durch die ausgezogene Linie dargestellt ist, als |f 


w, an w, angesetzt. Auf diese Weise wird bei r, die Durchbiegung zu Null, was ohne weiteres zu 

erkennen ist. Aber auch die Kriimmung verschwindet, weil sie bei der Drehung der abgeschnitte- 
/ . . 

nen Flache um 180° ihr Vorzeichen nicht andert und das Ansetzen von wj an w, einer Subtraktion 


gleichkommt. Der in Abb. 6 schraffierte Bereich gibt einen Schnitt langs der 7-Achse durch die | 


gesuchte Fliache wy. 


Die am Rande r, frei werdenden Randwerte von w; sind meist bedeutungslos und brauchen } 
in der Regel nicht mehr beseitigt zu werden; andernfalls ist bei einem drehbar gelagerten Rande | 


r, entsprechend wie bei r; zu verfahren. 


Ist dagegen der Rand rz einzuspannen, so sind fiir die dort durch den Schnitt frei gewordenen 


< ‘ We / ° ° 2 
Durchbiegungen und Verdrehungen neue Konstante eC, und JmC;, zu bestimmen. Die Durch 
biegungen w,, sind dabei der neuen Rechnung unverandert zugrunde zu legen, hingegen miissen 
die Verdrehungen mit umgekehrten Vorzeichen, als sie sich aus w, ergeben, versehen werden, weil 


fiir den Rand rz; nunmehr das Koordinatensystem (é’, 7’, w) gilt. Aus Abb. 6 ist dies ohne weiteres y 
zu ersehen. Die folgenden Beispiele geben einen Einblick in die Handhabung der beschriebenen | 


Methoden und der Zahlentafeln!. 


1e2 42 Ort dep 
| Singular itat re 


- 


aw, 
W473, an] / 


it eee 


Abb. 6. Ermittlung der Fliche w/; bei drehbarem Rand rz. Abb. 7. Lage der Koordinatensysteme beim ersten Beispiel. 


6. Beispiele. Zur Vereinfachung wird in den Beispielen die Querdehnungszahl vy = 0 gesetzt, 
ohne da®B dadurch der Rechnungsgang im Prinzip beeinfluBt wird. 


1) Im ersten Beispiel soll das EinfluBfeld eines Stiittzmomentes in der Mitte eines Randes 
einer allseits eingespannten quadratischen Platte ermittelt werden. Nach Pucher ist die Singulari- 
tat fir das Biegungsmoment am eingespannten Rand 

1 1”? 

aN &+72~ 

Sie verlauft symmetrisch zur 7-Achse (Abb. 7). Aus der Forderung, daB gemaB der endgiltigen 
Lagerung der Platte an den vier Ecken mit den Koordinaten (+1, 0) und (+1, —2) die Durchbie- 
gungen und Verdrehungen der Flache w, nach beiden Achsrichtungen verschwinden missen, 
bestimmt man nach (2) die Flache we und erhalt nach einfacher Rechnung 

n 2 13 2 4 
7 ae + 557 95 fe 95 f'n 25 n°. (16) 


Wp 


—n Nw, =—2 N(wp + We) 


a) Rand r;. Zunachst wird der Rand r,; mit Hilfe der Flache w, eingespannt. Zur Berechnung 
von w; ermittelt man daher fiir 7=—2 die Durchbiegungen und Verdrehungen in der 7-Richtung 
der F'lache (16), also die Funktionen f,; und @s3, und bezieht sie auf das Koordinatensystem 


' Im Heft 6 der Forschungshefte aus dem Gebiete des Stahlbaus, herausgegeben vom Deutschen Stahlbau- 


Verband Berlin, Berlin 1943, S. 189 bis 199 hat G. Worch in seinem ,,Beitrag zur Berechnung der allseitig ein- — 
gespannten Platte‘* das vom Verfasser entwickelte Verfahren, wie es in dem in FuBnote 3 von Seite 106 an- — 


gefiihrten Buch enthalten ist, mit geringen Abinderungen iibernommen und zur Berechnung zweier einfacher 
Belastungsfalle der quadratischen Platte benutzt. Da die Entwicklungen von Worch das Seitenverhaltnis ent- 


halten, besitzen sie nicht den Vorteil der besonders in der vorliegenden Arbeit dargelegten Tabulierungs- 


moglichkeit, 
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(€,.7,), was zur Folge hat, daB sich die Vorzeichen von sz umkehren: 


4 
aN fs, = TU IN 5525 Ss tee 0,96 + 0,16 &, 
dw 4&2 : 
Baer. N —— 0 a S1 cor 2 
mN gg = — 2. N (“2 : ie Gerbar 024 + 0.088. 


Da diese beiden Ausdriicke symmetrisch zur 7,-Achse verlaufen, sind die 'Tafeln des Symmetrie- 
falles maBgebend. : 
Es. werden nun f,; und @sz im Intervall 0 = €, = 1 in den Punkten 0, 0,1, 0,2...1 und mit 
+1 


41 
der Simpsonschen Regel die Integrale Zig= — | fe Re Da dé, und Zig = — =H) Psg Ite Dn dé, be- 
i il 


- rechnet (s. Tafel Aa). Dabei wird einmal h=1, das andere Mal h=2 gesetzt, also mit der zweiten 
Naherung gearbeitet. Man erhAlt Z,f=—0,024843, Z,,=0,163107, Z,f¢=0,004756, Z,,=—0,059226. 
Diese Werte werden in die Gleichungen fir die Konstanten %teC,, 3mC,, ReC,, mC, auf 
Tafel Ab unter I, zweite Naherung eingesetzt und ergeben Re C,= 0,043089, Sm C, = 0,041898, 
Re C, = —0,000463, Sm C, =0,000800. Mit diesen GréBen geht man in die Tafeln Ac 1 und Ac2 und 
gewinnt die Durchbiegungen der Flache w3. Die Tafeln Ad1 und Ad2 liefern die zugehérigen 


Verdrehungen, die aber nur langs der Rander r; und r, von Interesse sind. 


Es ist wichtig, die Genauigkeit der eben aufgestellten Rechnung zu tiberpriifen und festzu- 
stellen, wie weit die Randwerte f.; und q.3 der singularen Flache langs r, durch das Hinzufiigen 
von w; zum Verschwinden gebracht werden. Die folgende Zahlentafel 1 enthalt in der zweiten 


Zahlentafel 1 


E | 0 0,2 0,4 0,6 eee | 1 
—aNfeg +0,04000 +0,03650 | +0,02714 | +0,01503 | +0,00447 0 
w', —0,04081 —0,03626 | —0,02609 | —0,01598 | —0,00665 0 
—aNqs, | —0,24000 | —0,22699 | —0,19022 | —0,13545 | —0,06989 0 
(G,/I7), | —+0,23298 | +0,22885 | +0,20057 +0,12822 | -+0,03875 0 


und vierten Zeile die gegebenen GréBen —xzNf., baw. —aN@qsy, in der dritten und finften Zeile 
stehen die von w$ herriithrenden Randwerte w3, und (dw3/d7,),. Entsprechend (6) miissen die 


Gleichungen 


s 0 ws 
—2 Nfs +3, = 0, —aNoa + (a ),=9 


moglichst gut erfillt sein. Man sieht, daB dies besonders in der Umgebung der Symmetrieachse 
(€,=0) der Fall ist. Nach den Ecken (£,=0,8) werden die Abweichungen allerdings verhaltnismabig 
gréBer. Dies ist aber bedeutungslos, denn gegeniiber der gréBten Ordinate der endgiltigen 
EinfluBflache, die 1 ist, betragt bei €=0 der Fehler nur 0,00218, liegt also bei 2°/o). Nur wenig 
schlechter sind die Verhaltnisse bei den Fehlern der Randverdrehung. 

Die Flache w§ verlauft in der 7-Richtung nach einer gedampften Schwingung. Am Rande r, 
werden daher weder w$, noch (dw,/07,), zu Null, geniigen also nicht den Bedingungen der Ein- 
spannung. Da fir w3 das Koordinatensystem £,7, maBgebend ist, erhalt man die GréBtwerte von 
ws, und (dw3/d7,);, wenn man unter Beriicksichtigung der oben ermittelten Konstanten w; und 
(dw3/d7,) fir £,=0, 4, =—2 mit Hilfe, der Tafeln Acl und Ac2 bzw. Adl und Ad2 berechnet. 
Fuhrt man darauf wieder das Koordinatensystem (&,7) ein, so ist das Vorzeichen von (01w3/d7,)1 
umzukehren. Es ergibt sich dann 

» dws \ ze 
wi, = —0,00335, Sih — 0,01077. 
Diese GréfSen sind sehr klein und kénnten ohne Schaden fir die Genauigkeit der gesuchten 
HinfluBflache vernachlassigt werden. Zur Erlauterung der am Ende des Abschnittes 5 erwahnten 
Berechnung der Flache w’ soll aber auch diese von w} herriihrende Randstorung beseitigt werden. 

Die Tafeln Acl und Ac2 bzw. Ad1 und Ad2 liefern fir 7,=—2 die maBgebenden Randwerte 
von w langs r, zwischen 0 S &, = 1 nur in den Punkten &=0; 0,2; 0,4;...1. Um auch die Werte 
in den Zwischenpunkten = 0,1; 0,3; ...0,9 zu erhalten, interpoliert man am besten graphisch, 


y 
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was vollig ausreichend ist. Da die Randwerte w3, und (dw3/d), klein sind, kann. man sich bei 
der Ermittlung der Flache w3/ mit der ersten Naherung begniigen. Man berechnet also ebenso 
wie oben die Werte 


Z,7 = 0,002033, Live ——0,006638 | 
mit Hilfe der Tafel Aa und darauf die Konstanten 
Re C, =—0,001979, ym C, =—0,002284 


mittels der Tafel AbI, erste Naherung. Bei €=0 erhalt man dann w%,=0,00288 gegeniiber 
wi, =—0,00335 und (dw3’/d7),=—0,01160 gegentiber (0 w3/d 4), = 0,01077. Die Abweichungen 
sind also in Anbetracht der Kleinheit der gegebenen Werte ohne Bedeutung. 


b) Rand r,. Wir wenden uns nun den Randwerten langs r, zu, die durch wy erzeugt werden. 
Dort ist wo, = 0, wahrend gewisse Verdrehungen (dw,/dy), vorhanden sind. Diese werden wiederum 
zwischen €=0 und &=1 in Abstanden von 0,1 aus 


nN (5%) =a Ng = — 0,08 (1-8) aN 
n=0 


oi) 
berechnet. Mit Hilfe der Tafel Aa erhalt man 
Lif mZof=s Z19 =—9,056288 , Zoey = 9,023869 


und weiter mit Tafel Ab 
Re C, =—0,012322, mC, =—0,009394, Re C, = 0,000465, Sm C, =—0,000048. 


Die Zahlentafel 2 gibt eine Ubersicht tiber die Fehler. Das Eingehen auf die Restwerte dieser 


Flache w, am Rande r, eriibrigt sich, da sie unbedeutend sind. 


Zahlentafel 2 


E 0 0,2 0,4 0,6 0,8 1 
—naNfry 0 0 0 0 0 0 
wy +0,00045 —0,00014 —0,00055 +0,00053 +0,00107 0 
—xnN@gs, |° —0,08000 —0,07680 —0,06720 —0,05120 —0,02880 0 
(dw3/d7), +0,07580 +.0,07799 +0,07304 +0,04675 +0,01283 0 


c) Rander r, und ry. Es bleibt nun noch die Aufgabe, die Rander r, und r, einzuspannen. 
Wegen der Symmetrie der Flache wy zur Achse 7 werden beide gleich behandelt. Man ermittelt 
als erstes die Durchbiegungen und Verdrehungen in der €-Richtung der Flache wy am Rande ry, 
(&=—1) aus (16). Es ergibt sich 


_aNf,=—nNw,__, =) a 0,52? — 0,163, 
oe sys pe NIPENE: 
mae Ds =-aN (5 Poe (+1) Spas 


Beide Funktionen werden im Intervall 0 <= y < —2 zahlenmaBig in Abstanden von 0,1 berechnet. 
Unter Zugrundelegung des Koordinatensystems (£,72) zerlegt man dann die unsymmetrischen 
Funktionen fy und gy, in die Symmetrieteile fs. und gsy und die Antimetrieteile fag und Qap. Es ist 
wieder zu beachten, daB im System (£5,179) die Vorzeichen der Verdrehungen gs. und az gegen- 
tuber gy, umzukehren sind. Die Berechnung der Integrale erfolgt fiir den Symmetrieteil nach 
Tafel Aa: 


LZ, f= —9,087797, Z,9 = 0,223836 , 237 = 0,01 1950, Z og =—9,064488 . 


Die zugehorigen Konstanten sind nach Tafel Ab 
Re C, = 0,074661, 3m C, = 0,089200, ReC, = 0,000644, Sm C, = 0,002095. 
Die entsprechenden Werte des Antimetrieteils erhalt man nach Tafel Ba und Bb: 


Zip=  9,011603, 1g =—0,023834, Z3¢ =—0,003544,.  Z3 4 =—0,002850. - 
Rte C*¥ =—0,007078, YmCT =—0,005346, ReCZ =—0,000979, GmCz =—0,000831. 


Bi Mt eee ee rn ae ; pee 
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Zahlentafel 3 


§ 


0 


0,2 


0,4 


0,6 


0,8 


—aNfry +0,14000 +0,12860 +0,09703 +0,05412 +0,01571 

ws, —0,14073 —0,12843 —0,09597 —0,05496 —0,01830 
—2N sy —0,34000 —0,31982 —0,25834 —0,16090 —0,05721 
(dw$/d 775)» +0,34115 +0,31922 +0,25787 +0,16253 +0,05354 


—2Nfos +0,01076 +0,01502 +0,01086 +0,00324 
w%, 0 —0,01028 —0,01576 —0,01055 —0,00215 
—2 Nap 0 —0,02808 —0,03493 —0,01291 +0,01524 | 

(dw2/0 72)> 0 +0,03605 +0,02301 +0,01503 +0,01724 


— a) — a) fal) — ee 


Die Zahlentafel 3 gibt einen Einblick in die GréBe der Fehler. Wahrend wi, (0.w3/d 7). und w’, 
in recht befriedigender Weise die entsprechenden gegebenen Randwerte fs», Qs. und fa, zum Ver- 
schwinden bringen, fallen die starken Unterschiede bei —aNq@az und (dw$/dn.), auf. Das liegt 
daran, daB die Funktion —NGay zwischen &,=0 und &=1 baw. £.=0 und £,=—1 ihr Vorzeichen 
wechselt (Abb. 8). Diese starke Schwankung kann aber durch die zweite Naherung des Anti- 
metriefalles noch nicht gut genug angeglichen werden. Die dritte Naherung, deren Rechnung 
hier im einzelnen nicht wiedergegeben ist, ergibt bessere Werte, die in Abb. 8 durch die 
gestrichelte Linie gekennzeichnet sind. 


a 
(oe) 
On2 /2 


Abb. 8. Graphische Darstellung der antimetrischen Verdrehungen 
der Zahlentafel 3. 


Abb, 9. Graphische Darstellung der Verdrehungen 
der Zahlentafel 4. 


Zahlentafel 4 


e ps Sa 0 Od nO ern Oh sg 10,0" OAT Oe he oe ea et 
—aNfz 0 | +0,01247 | +.0,04326 | +0,08200 | +.0,11784 | +0,14000 | +0,13936 | +0,11207 | +-0,06497 | +0,01894| 0 
Wo 0  |—0,01615 |—0,04441 |—0,08021 |—0,11815 |—0,14073 |—0,13871 |—0,11173 |—-0,06551 |—0,02045| 0 

if 
—aN 0  |—0,07245 |\—0,14799 |—0,22341 |—0,29174 |—0,34000 |—0,34791 |—-0,29327 |—0,17381 |—0,04196| 6 
(w2/d72)o} 0 | +0,03630 | +0,14750 | + 0,23486 | +0,28317 | + 0,34115 | + 0,35527 | + 0,28088 | +0,17756 | +0,07078| 0 


Die eben besprochenen grofen Fehler von (dw$/d72). werden allerdings sehr stark abge- 
schwacht, wenn man Symmetrie- und Antimetrieteile von f, und qm, addiert, weil die Antimetrie- 
teile nur einen Bruchteil der Symmetrieteile ausmachen (Zahlentafel 4). Entscheidend fiir die 
Beurteilung der Brauchbarkeit der Rechnung sind aber allein die Unterschiede zwischen den 
gegebenen unsymmetrischen Funktionen —aNfy baw. —aNqg, und ihren Annaherungen Wye 
bzw. (02,/072)>. Diese betragen tiberall mit Ausnahme der Nachbarschaft der Ecken nur wenige 
Prozent (Abb. 9). Die Differenzen an den Ecken haben auf die Form der endgiltigen Flache 
keinen wesentlichen Einflu8, wie bereits erértert wurde. 

Die Restwerte der Flachen w$ und w$ am gegeniberliegenden Rande r, sind klein und kénnen 
vernachlassigt werden. Beabsichtigt man aber, auch sie zu beseitigen, so ist mit ihnen ebenso wie 
mit den Restwerten der Flache w§ am Rande r, zu verfahren. 
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Die einzelnen, eben berechneten Teilflachen werden mit wy zusammengesetzt; wobei die 
demselben Rand zugeordneten Konstanten der Symmetrie- baw. Antimetrieteile vor ihrer 
Multiplikation mit den Flachen der Tafeln 
A und B addiert werden kénnen, wodurch 
die Zusammensetzung tbersichtlicher wird. 
Die Abb. 10 gibt ein Bild der gesuchten — 
EinfluBflache in Schichtenlinien?. 


2) Das zweite Beispiel behandelt eine 
quadratische, gleichmaBig belastete Platte, 
bei der zwei Nachbarrander eingespannt und 
die beiden anderen drehbar gelagert sind. 


— 
71 
Abb. 10, EinfluBfeld fair das Stiitzmoment im Punkte €=0, 7=0 Abb. 11. Lage der Koordinatensysteme 
einer allseits eingespannten quadratischen Platte, 8 a-fach. beim zweiten Beispiel. 


Man erganzt das in Abb. 11 mit ausgezogenen Linien dargestellte quadratische Plattenfeld 
nach der negativen Seite der &-Achse zu einer Platte mit den Abmessungen 2-1 und belastet 
sie antimetrisch zur Achse 7 mit der konstanten Last p. Auf diese Weise wird erreicht, daB 
langs £=0 Durchbiegungen und Momente der zugehérigen Biegeflache verschwinden, und daB 
dort die Verhaltnisse wie an einem drehbar gelagerten Rande herrschen. Die Gleichungen fir 
Ww, sind 


24.N 24.N : 

oe om pete to) — eh OE EST aaa) | 

24.N 24 N ‘ ete te, (17) 
pa Wo = pat (wp 4 We) = —§ +0,5€—1,5€ (f= 0)5 


Es ist bei dieser Flache 0 w,/0 7 = 0. Mit Hilfe der Tafeln Ba und Bb erhalt man nach Berechnung 
der Funktionen (17) in den Zehntelpunkten zwischen €=— | und &=1 die Integrale 


Zi} = 0,102642, Z3¢ =—0,040437, Lig=Lig=9 
und die Konstanten 


Re C¥ =—0,035174, QmCf =—0,036931, Re C¥ =—0,000798, QJmCz =—0,001998. 


Zahlentafel 5 


é 0 0,4 
ta 0 0,08960 0,12960 0,10560 0,04160 0 
wt, 0 —0,08937 —0,13005 —0,10516 —0,04205 0 
Pat 0 0 0 0 0 0 
(dw%/d); 0 +-0,00136 —0,00202 +0,00257 —0,00376 0 


Zahlentafel 5 zeigt die recht geringen Unterschiede der Randwerte der so berechneten Flache wf 
gegentiber den gegebenen der Flache wy. 

Am Rande rz (y=—1) ist wf durchzuschneiden. Da r, drehbar gelagert werden soll, eriibrigt 
es sich gemaB den Erérterungen im Abschnitt 5, fiir die durch den Schnitt frei gewordenen 
Randwerte neue Konstanten zu berechnen, wie es im ersten Beispiel bei den eingespannten 


1 Die Abb. 10 ist dem in FuGnote 2 von Seite 106 erwahnten Aufsatz von A, Pucher, Ing.-Arch. 12 (1941), 


S. 91, entnommen. 
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Randern erforderlich war. Es wird vielmehr der mit den gleichen Konstanten wie oben multi- 
plizierte Teil der Antimetrieflachen, der zwischen 7=—1 und 4 = —2 liegt (s. die Tafeln B), von 
dem zwischen 7 =0 und 7 = —1 liegenden entsprechenden Flachen abgezogen (s. Abb. 6), und zwar 
so, daB 7=—2 mit 7=0 zur Deckung kommt. Danach ist z. B. die aus wt und w%’ herrihrende 


Durchbiegung im Punkte €=0,2, 4 =0,4 . 


—-F (wy + wi) = —0,035174 (—0,01013 — 0,00680) — 0,036931 (0,61270 + 0,00050) — 
— 0,000798 (0,14454 — 0,00008) — 0,001998 (0,41346 + 0,00007) . 


Die Zahlen innerhalb der Klammern sind den Tafeln Bel und Be2 entnommen. 

Die Randwerte, die von w*’ am Rande r, ubrig bleiben, sind so klein, daB sie nicht mehr 
beseitigt zu werden brauchen. 

Nunmehr ist der Rand r,; der Flache wo unter Zugrundelegung des Koordinatensystems 
(€,,7,) drehbar aufzulagern. Man verfahrt im Prinzip genau wie beim Rande r,. Die Durch- 
biegungen der Flache wy sind bei r, die gleichen wie bei r,, da wo eine Zylinderflache ist, ferner 
ist auch (0? w 9/073) = 0. Man erhalt 

if = —0,102642., 3¢ = 0,040437, Lig Liew 220: 


Die Vorzeichenanderung dieser Integrale gegentiber den entsprechenden am Rande r, ist auf die 
gedrehte Lage des Koordinatensystems (¢,,7,) zuriickzufiihren. Nach Tafel Bb IL, zweite Nahe- 
rung ergibt sich 

Re CT = 0,010358, SmCy =0,024620, ReC = —0,000059, 3mCz = 0,000444. 


Beim Gebrauch aller Tafeln B des Antimetrieteils ist in bezug auf die Vorzeichen stets darauf 
zu achten, wie das jeweils giltige Koordinatensystem zu den Antimetrieflachen wy und w* sowie 
den Ableitungen von w* angeordnet ist. 


Zahlentafel 6 gibt einen Uberblick wher die Genauigkeit der Rechnung. 
Zahlentafel 6 


é | 0 | 0,2 0.4 0.6 0,8 1 

a 0 0,08960 0,12960 0,10560 0.04160 0 

w', 0 —0,08927 —0,13019 | —0,10538 | —0,04061 0 

Was 0 | 0 Seah 0 0 0 
(8229/6 73), 0 | —0,00642 | +0,00884 | +0,00164 | —0,03650 0 


Am Rande r, hinterlaBt der dort zu fiihrende Schnitt 
durch wf freie Durchbiegungen und Verdrehungen, die aus 
den Tafeln Bel und Be2 bzw. Bdl und Bd2 berechnet 
werden. Da r, einzuspannen ist, sind neue Integrale und 
Konstanten zu ermitteln. Diese sind 


Z¥;= —0,003361, t= 0,015082, 
Z3,—0,000991, Z5 , = —9,004966 


Re CT = 0,002971, Sm CT = 0,002038, 
Re Cl = 0,000228, Sm CZ = 0,000264. 


Die zugehérige Flache ws kann am Rande rz ebenso ,,ge- 

spiegelt’* werden, wie es bereits bei wi geschehen ist, so dal 
. . oe Wi . 

man ohne neue Rechnung eine weitere Flache wi gewinnt. 

Es sollen:- nun noch die Stiitzmomente und die Feld- 

A i i n arallel Zu den Ko- Abb. 12, Biegungsmomente einer gleichmaBig 

gegmLcnte langs eae Diagonale P belasteten quadratischen Platte at zwel ie 


ordinatenachsen bestimmt werden. Die Gleichungen fiir die — nachbarten eingespannten und zwei drehbaren 
j j Randern. 
Biegungsmomente sind 


und 


> 


WAY he 
/ thee pa 


N / 0? w 02 w N / 0? w 0? w 
( aoe ae mM, = 5 =e eran 
Wes a7? a\on dé 
Wie bereits eingangs erwahnt, wird » = 0 gesetzt. Die Momente sind symmetrisch zu der in Abb. 12 
eingezeichneten Diagonale. Es ist vorteilhafter, die Einspannmomente am Rande r, zu berechnen, 
weil hier genauere Werte als bei r, zu erwarten sind, denn bei r, erfiillen samtliche Teilflachen die 


Randbedingungen identisch. 


Me = ae 
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Die einzelnen Anteile, welche die verschiedene 
und die Endergebnisse sind in Z 


Koepcke: 


die dem gleichen Rand zugeordnet sind, 


liegende Beispiel die be 


und 


iden Gruppen 


ahlentafel 7 und 8 zusammengeste 


Re C* = —0,035174 + 0,002971 = 20,032203, 


3m C¥ = —0,036931 + 0,002038 = —0,034893, 
Re CY = —0,000798 + 0,000228 = —0,000570, 


3m CE = —0,001998 + 0,000264 = —0,001734 


ReC* =+ 0,010358, 
Re C} = —0,000059, 


gm C¥ = 0,024620, 
Sm CX = 0,000444 


Bei der numerischen Berechnung der Momente ist wiederum den Vorzeichen der 


Anteile besondere Aufmerksamkeit zu schenken. 


Zahlentafel 7 (Stiitzmomente) 


EinfluBflachen und innere Krafte umfangsgelagerter Rechteckplatten. 


zugehorig zu ry 


zugehérig zu r3- 


he YAS AS 7%) <<. 
NAAR Nie eee gy 
2g lan rp pate 


Soh 
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n Biegeflachen zu den Momenten beisteuern, 
lit. Man kann alle Konstanten, 
addieren und gewinnt auf diese Weise fiir das vor- 


einzelnen 


—0,2 —0,4 —0,6 —0,8 a 
aus Wo —0,12500 —0,12500 | _—0,12500 —0,12500 —0,12500 | —0,12500 
aus w?tw*’ | +0,13371 +0,09198 +0,05157 +0,02566 +0,01164 | +0,00481 
aus wi’ tw} +0,00006 +0,00004 —0,00009 —0,00052 —0,00176 | —0,00481 
aus w$ +0,00192 +0,00544 +0,01352 +0,03051 +0,06286 | +0,11643 | 
més 0 —0,02754 —0,06000 —0,06935 —0,05226 | 0 | 
(+0,01069) (—0,00857) 


$e 


ate 3 


a 


| 


| 
| 
| 
| 


Zahlentafel 8 (Feldmomente) 

E 1 0,8 0.6 0,4 0,2 0 

M7 ede [ee a) —0,4 ie iG eS et 
aus Wy —0,12500 —0,02000 | +0,04500 -+0,07000 +0,05500 0 
aus wt +w$’ +0,13371 +0,01676 —0,02147 —0,02056 —0,00817 0 
aus wt’ +-we’” | +-0,00006 —0,00008 +0,00007 +0,00096 | +0,00189 | 0 
aus w% +0,00192 +0,00295 —0,00644 —0,02530 —0,03328 | 0 

mef 0 —0,00037 +0,01716 +0,02510 +0,01544 ¢ 

(+.0,01069) 


Die Zahlentafel 7 enthalt die Stiitzmomente m; langs des Randes ry, in der Zahlentafel 8 sind 
die Momente m; langs der Diagonale ermittelt und in Abb. 12 aufgetragen. Es fallt in den beiden 
Tafeln auf, da® die Endsummen der Biegungsmomente in den Ecken nicht Null werden, was 
aber notwendig ist, weil die vier Rander eine Ebene bilden. Die Abweichung hat ihre Ursache | 
darin, da®B langs r, und r; noch kleine Restdurchbiegungen verblieben sind, die wellenformig 
verlaufen. Diese Wellen erzeugen Momente, deren Einflu8 allerdings im Plattenfeld sehr bald 
verschwindet. Man kann daher sicher sein, da® die errechneten Biegungsmomente schon in 
geringer Entfernung von r, und r; ausreichend genau sind. Die sich aus der Rechnung ergebenden, 
von Null abweichenden Werte in den Ecken werden deshalb unterdrickt und an ihre Stelle die 
richtigen Werte gesetzt. . 


(Eingegangen am 23. April 1949.) 


_ Anschrift des Verfassers: 
Dr.-Ing. Werner Koepcke, (1) Berlin-Charlottenburg, Technische Universitat, Lehrstuhl fiir Stahlbetonbau. 
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Nachtrag 
zu meinem Aufsatz in Band XVII, S..187 des Ingenieur-Archivs: 


, Die energetische Berechnung des induzierten Widerstandes“. 


Von W. Kaufmann. 


Herr E. Truckenbrodt! hat mich freundlicherweise darauf aufmerksam gemacht, das die in 
Abb. 3 meiner obigen Arbeit aufgetragene Funktion y?/n=f(n) bei elliptischer Fligel-Auftriebs- 
verteilung fir 7 + 0 nicht gegen oo geht, sondern einen endlichen Grenzwert besitzt. 

Kine besondere Entwicklung in der Umgebung der Stelle 7=0 ist also in diesem Falle nicht 
notwendig. Der Ausdruck fiir den induzierten Widerstand lautet dann unter Zusammenfassung 


der beiden Anteile aus Gleichung (8) und (15) einfacher 


n=1 

> Ole v 

Wi, = &" Cas [ San), (1) 
n=90 


_wobei das Integral durch graphische Integration tber den ganzen Bereich 0S 7 = 1 bestimmt 
werden kann. Obwohl sich dadurch an dem von mir mitgeteilten Zahlenwert praktisch nichts 
andert — auch Herr Truckenbrodt berechnet nur eine einprozentige Abweichung gegentiber dem 
Prandtlschen Wert —, méchte ich nachstehend doch die von Herrn Truckenbrodt durchgefihrte 
Grenzwertbetrachtung in etwas erweiterter Form wiedergeben, da sie mir fiir die Anwendung des 
Verfahrens von Bedeutung erscheint. : 

Unter Beachtung der Gleichungen (2) und (18) meiner Arbeit und nach Einfithrung der Ab- 


= 


kurzung c=(«,—x)/a kann Gleichung (17) wie folgt geschrieben werden: 


(a) ae) (VIR ee Les) 


oder wegen (3) 


iy a eee | 
7 SS pS 
Sg (V14+2e?—1+4c%)  (0XeS 1). 
Durch Differentiation nach ¢ folgt daraus l 
— + 1) Cc 
dy 1 V1+2¢ 
Core Wii roe Sine 


und durch Grenzwertbildung 


di a 
| ae ee ay V3 


Nun wird wegen c=(x,—x)/a und unter Beachtung von Gleichung (16) und (18) der obigen Arbeit 


4 4 
py [xy 2 Foe Te 2P a 
Vevey hon (=) gal ee Ata Sighs 20K 
lim ( = - | 4, be c ie x—>b/2 E=6/2 
q if Tad 
3 3 
iat SERVE the fet eV gash | rey AUR 
x—>b/2 ede x—>b/2 pedi ities 
Jf vedas Sf vedé 
ee E=b/2 
Da aber ae 
d(yx*) 2 1x de Lop 
et te: 
£=b/2 
so folgt ed 
ey 3 Ve ee a eta E aa 
lim (7) seg! 4 dx |x=b/2 es Saale sig Meares 


wobei yx=I'x/I, die dimensionslose Fliigelzirkulation darstellt. Fihrt man noch die dimensions- 
losen Werte x=2 x/b und a=2a/b ein, so erhalt man schlieBlich den von Herrn Truckenbrodt 


1 Anschrift: Dr.-Ing. E. Truckenbrodt, (20b) Braunschweig, Andreeplatz 5. 


140 Kaufmann: Nachtrag. 
2 
gefundenen Grenzwert dp | 
ey SW ee —[d ve) | 40 | 
me ee Fb 13208 dx |(z=1) \ henstante. Zirkulations- | 


pecs verteilung . 
Im Falle elliptisch verteilter Flagelzirkulation 


Haley | x2 wird unter Beachtung von (2a) 


45 


lim (= = 2,19G = 2,79 = 2,19. 

n—9 1 oi 
Herr Truckenbrodt hat dariiber hinaus den induzierten 
Widerstand fir elliptische, parabolische und dreiecks- 
formige Verteilung der Flagelzirkulation nach dem 
vorliegenden Verfahren berechnet und mit den ent- 
sprechenden Werten nach der Prandtlschen Theorie 


verglichen. Die Ergebnisse seiner Rechnung sind in 


(pacaalieye 
f Zirkulationsverteiling 


Grerecktormge 
ZirkulotionsverteHung 


der Tabelle und in nebenstehender Abbildung zu 5 22 28 26 28 0 
sammengestellt. Aus der Tabelle ist die gute Uberein- ‘ 
stimmung beider Theorien erkennbar. Abb. 1. Zuordnung (n) fair verschiedene Zirkulationsverteilungen, 
Zirkulationsverteilung | | Induz. Widerstand W; | 
; - d(yx') 
x S| Prandtl K 
Vs a i= ee ran CAME ve Wi 
. er. abl = = = 
Noth 72 % be wo ol Ww. ol 
elliptisch 1—x Z 2 hee 1,01 mone 
. 3 ene 2 : 2 
parabolisch [eon — 0 ged 1,01 ely 
3 cs I 
dreiecksformig 1—% an 0 In patios 1,01 In 9, of o 
2 us It 
Bei beliebig vorgegebener, zwischen Ellipse und umbeschriebenem Rechteck liegender Zirku- 


lationsverteilung kann man, wie ich an anderer Stelle gezeigt habe’, die Flagelzirkulation an- 
genahert durch den Ausdruck darstellen Sea: ‘e 
r.=T,|1— (=) | 
wo y>2. In diesem Falle geht [y?/7],-0 > ©, 8° da® man zweckmaBig Gleichung (20) meiner 
obigen Arbeit benutzen wird. 

Wirde man dagegen, worauf Herr Truckenbrodt hinweist, die Zirkulationsverteilung durch 
den Betzschen Ansatz rate 


T= 1; - Gey > An (3). 


darstellen, so wiirde sich [y2/y],—0 stets endlich ergeben, solange > A, + © wird ; In diesem Falle 
wirde die obige Gleichung (1) anzuwenden sein, da jetzt die graphische, Integration iiber den 


ganzen Bereich 0 <7 1 keine Schwierigkeiten bereitet. 

Die kritischen Bemerkungen des Herrn Truckenbrodt und die von ihm durchgerechneten Bei- 
spiele stellen eine wesentliche Ergainzung meiner Arbeit dar, fiir die ich ihm aufrichtigen Dank sage. 

Zum Schlu8 méchte ich noch einen Druckfehler in meiner Arbeit richtigstellen. Auf 5. 192 | 
Zeile 3 mu es, wie man leicht feststellt, richtig heiBen ; 

ee) i 
a] 
1 W. Kaufmann, Sitzungsber. Bayer. Akad. Wiss. Mathem.-Naturwiss. Abt. 1946, 5. 119. 
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(Eingegangen am 27. September 1949.) 
Anschrift des Verfassers: Professor Dr.-Ing. Walter Kaufmann, (13b) Miinchen, Technische Hochschule. 
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‘der Supraleitung 


Von 
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Zweite Auflage 
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Finleitung. 

Diese Schrift strebt die Klarung der Anschauungen iiber Supraleitung an, und zwar 
durch Ausdehnung der Maxwellschen Elektrodynamik auf die Supraleiter nach Ideen, 
welche Fritz London und Heinz London 1935 und spater angegeben haben. Sie soll genau 
soweit gehen, wie diese Erweiterung die Tatsachen deutet. Wie groB der von ihr umfaBte 
Bereich ist, wird der Leser sehen. Aber dariiber hinaus gibt es eine gewisse Zabl abweichen- 
der, dem Klaren Verstandnis bisher verschlossener Tatsachen. Wenn wir sie ignorieren, 
so berufen wir uns dabei auf Clark Maxwell, welcher die seit Newton wohlbekannte Dis- 
persion des Brechungsindex nicht in seine elektromagnetische Lichttheorie einbaute, 
auf Heinrich Hertz, der in seiner Elektrodynamik der bewegten Kérper den Fresnelschen. 
Mitfihrungskoeffizienten nicht beriicksichtigte, und auf jene Vater der Thermodynamik 
und Gastheorie, weléhe von der einfachen Zustandsgleichung pV = RT. ausgingen, ob- 
wohl sie wuBten, daB diese eigentlich fiir kéin Gas gilt. Die Theorie sucht sich eben der 
Wirklichkeit mittels sukzessiver Approximation zu nahern. Zudem haben wir bei mancher 
unerklirlichen Tatsache den Verdacht einer gewissen Unvollkommenheit des Versuchs- 
materials. Der ideale feste Kérper ist ein Einkristall; die meisten Supraleitungsversuche 
aber sind mit polykristallinem Material angeste'lt. Und selbst der beste Einkristall, den 
je ein Mensch in der Hand hatte, ist, wenn er ein paar Zentimeter grof ist, ein ,,Mosaik- 
Kristall, ein Gefiige aus vielen kleinen Kristalliten, an deren Grenzen das Raumgitter 
sich nicht stetig fortsetzt. Fiir die fast immer auftretenden Hysterese-Exscheinungen beim 
Ubergang vom Supra- zum Normalleiter und umgekehrt kommt hinzu, daB man auch 
sonst bei tiefen Temperaturen und insbesondere bei festen Kérpern haufig Verzégerungen 
findet, welche das theoretische Gleichgewicht zweier Phasen unter Umstanden vollig ver- 
decken. SchlieBlich-ist manche der nicht gedeuteten Beebachtungen nur selten, vielleicht 
nur einmal gemacht, so da8 Nachpriifung unter variierten Bedingungen notwendig er- 
scheint, bevor man dariiber nachdenkt. Die Supraleitung wurde 1911 entdeckt; aber trotz- 
dem ist ihre Erforschung vielleicht tiber das Anfangsstadium noch nicht hinaus, weil das 
Experimentieren bei so tiefen Temperaturen schwierig und nur an wenigen Stellen méglich 
ist. Bewahrt sich dieim folgenden vertretene Theorie, so leistet sie fiir den Supraleiter 
dasselbe wie die urspriingliche Maxwellsche Theorie fiir den Normal- und den Nichtleiter. 
Mehr will sie nicht. 


Die atomare Theorie der Supraleitung scheint uns‘noch in keiner ihrer mannigfachen 
Formen so weit geklart zu sein, daB wir sie in dieses Buch aufnehmen méchten. Immerhin 
steht wohl fest, da®B derim Prinzip ,,ewige™ Dauerstroi, den wir in einem supraleitenden 
Ringe herzustellen vermégen, das makroskopische Analogon zu jenen atomaren Dauer- 
stromen darstellt, welche das Bohrsche Atommodell seit jeher im Innern der meisten Atome | 
und Molekiile annimmt. In beiden Fallen liegen typische Quanteneffekte vor, die der 
alteren Physik unverstandlich bleiben. Gerade in dicsem unbestreitbaren experimentellen 
Nachweis solcher Quantenstréme scheint mir das Hauptinteresse der Physik an der 
Erforschung der Supraleitung zu liegen. Dieses Buch aber soll zeigen, wie sich diese 
Stromung in die Maxwellsche Theorie einordnen, 1aft. 
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In diesem Buch werden Theorie und Praxis der Destillier- und Rektifiziertechnik knapp, 

klar und dem neuesten Stand der Wissenschaft entsprechend dargestellt. Es ist fir den 

lernenden wie fiir den in der Praxis tatigen Ingenieur geschrieben, und seine besondere 

Starke liegt. in der gliicklichen Vereinigung strenger Wissenschaftlichkeit mit ingenieur- 

miBiger Behandlung des gesamten Gebietes, das in jiingster Zeit an praktischer Bes ; 
deutung noch zugenommen hat. : 


Die erste Auflage (1940) war schon kurz nach Redman vergriffen. Inzwischen ihed 
weitere Forschungen, nicht zuletzt des Verfassers selbst und seiner Karlsruher Schule, 
eine Reihe von Liicken in der Erkenntnis geschlossen. Die zweite Auflage, die wegen 
zeitbedingter Schwierigkeiten erst jetzt veranstaltet werden konnte, unterscheidet sich 
daher nicht unwesentlich von der ersten. So ist die Theorie des gekoppelten Warme- 
und Stoffaustauschs ausgebaut und vereinfacht, die u. a. fiir die Berechnung von Destillier- 
und Rektifizierapparaten und von Fullk6rpersiulen von groBer Bedeutung ist. Ferner 
haben die Vielstoffgemische gegeniiber der ersten Auflage eine wesentlich eingehendere 
Behandlung erfahren. Betrachtlich erweitert wurde der Abschnitt iiber Fiillkérpertechnik, 
und ganz neu aufgenommen ist eine Betrachtung tiber Sonderbauarten von Rektifizierbéden, 


So wird das Buch auch in der neuen Auflage seinen Rang als fiihrendes Werk auf diesem 
Fachgebiet behaupten. Sein Wert fiir die Praxis wird noch dadurch gesteigert, daB8 der 
Umfang des Anhangs, in welchem in Tabellen und Kurven Gleichgewichtswerte fiir — 
Zwei- und Dreistoffgemische und damit zuasammenhingende physikalische Werte wieder- 
gegeben werden, gegenuber der ersten Auflage auf ein Vielfaches angewachsen ist. 


Tnhaltsubersiehe: : \ 


Allgemeines. — Theoretische Grundlagen. — Flissigkeitstrennung dark einmalige De. 

stillation. — Die Rektifiziersaule. - Der stetig arbeitende Rektifizierapparat mit Aus- = __ i 
tauschbéden. — Die Rektifikationsvorgange im Warmeinhalt-Konzentrationsbild. — Tren- 
nung von Gemischen mit mehr als zwei Bestandteilen. — Bestimmung der Abmessungen 
der Rektifiziersiule mit Austauschbéden. Wirkung von Rektifizierbéden. — Rektifikation 
in Fiillkérpersaulen. — Flissigkeitsinhalt einer Rektifiziersiule. — Ausfihrung von Einzel- 
und Zubehérteilen. - Anhang: Gleichgewichtswerte fir Zwei- und Dreistoffgemische. — 
Namenverzeichnis. — Sachverzeichnis. 
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